Black-hole radiation in Bose-Einstein condensates 
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We study the phonon fluxes emitted when the condensate velocity crosses the speed of sound, i.e., 
in backgrounds which are analogous to that of a black hole. We focus on elongated one dimensional 
condensates and on stationary flows. Our theoretical analysis and numerical results are based on 
the Bogoliubov-de Gennes equation without further approximation. The spectral properties of the 
fluxes and of the long distance density-density correlations are obtained, both with and without an 
initial temperature. In realistic conditions, we show that the condensate temperature dominates 
the fluxes and thus hides the presence of the spontaneous emission (the Ifawking effect). We also 
explain why the temperature amplifies the long distance correlations which are intrinsic to this effect. 
This confirms that the correlation pattern offers a neat signature of the Hawking effect. Optimal 
conditions for observing the pattern are discussed, as well as correlation patterns associated with 
scattering of classical waves. Flows associated with white holes are also considered. 

PACS numbers: 03.75.Kk, 04.62.+V, 04.70.Dy 



I. INTRODUCTION 



The analogy between sound propagation in nonhomo- 
geneous media and light propagation in curved space- 
times has opened the possibility to detect the analogue 
of black hole radiation in the lab 1]. Indeed, when 
sound waves propagate in a flowing medium whose ve- 
locity crosses at some point the speed of sound, they ex- 
perience the analogue of an event horizon. If the phonon 
state is stationary and regular across this sonic horizon, 
one expects to obtain a thermal flux of phonons with a 
temperature fc^T = hK/2TT, where k is the gradient of 
the flow velocity evaluated at the sonic horizon. Since 
the analogy works perfectly in the hydrodynamical limit, 
the above result should be valid at least when k, is much 
smaller than the critical wave-vector characterizing the 
dispersion 

Following the original work of Unruh, various setups 
were p roposed, s ee B for a review. In Refs. 0, H, [1, 
[13, [H [13, [H, El [Hi the particular case of sound waves 
in dilute Bose-Einstein condensates (BEC) was consid- 
ered. These condensates have nice properties both from 
an experimental and a theoretical point of view. From 
the first, progress in the manipulation and control of 
their physical properties is rapid, and from the second, 
the equations for the condensate and the phonons are 
well understood, as well as the approximations they in- 
volve [3. 

In this work, we first aim to derive the properties of the 
phonon fluxes without making use of the gravitational 
analogy. In fact we also aim to determine the validity 
range of the analogy. To achieve these goals, our ana- 
lytical and numerical analysis are both based on the (ex- 
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act) Bogoliubov-de Gennes (BdG) equation. More specif- 
ically, we consider one dimensional stationary flows which 
contain one sonic horizon (black or white) surrounded 
by two infinite asymptotically homogeneous regions. In 
this case, at fixed conserved frequency \uj\, three types of 
asymptotic phonons exist, and a complete description of 
their scattering is given in terms of a 3 x 3 Bogoliubov 
transformation [l7j . 

Our second aim is to provide quantitative estimates of 
the spectral properties of the fluxes, in order to guide 
or explain experimental projects. To this end, we have 
performed a systematic numerical analysis. 

Our third aim is to understand the links between lo- 
cal and nonlocal observables. In Gaussian states (e.g. 
vacuum and thermal states) , the physical predictions en- 
coded in a Bogoliubov transformation are all contained in 
two groups of expectation values: first occupation num- 
bers (in the present case there are three of them, one 
for each type of outgoing phonons: = (a°"*^a°"*), 
i — 1,2,3), and second, interference terms, such as 
(a°"* a™j) with i ^ 7, which determine the long distance 
correlations [ll[ll,[2i. In this study, we were motivated 
by the recent work |14| where a distinct pattern of den- 
sity correlations was "numerically observed" when a sonic 
horizon forms. In what follows, we shall explain both the 
properties of this pattern, and also why, when taking into 
account the condensate temperature, the initial distribu- 
tions of phonons tend to hide the black hole radiation 
in the occupation numbers whereas, at the same time, 
they reinforce the correlation pattern without affecting 
its spatial properties. This confirms that the correlation 
pattern offers a neat signature of the Hawking effect [l3| . 

In Sec.[lll we derive the mode equation for the phonon 
field in one dimensional stationary condensates. We ob- 
tain an explicit fourth order equation which is valid for all 
nonhomogeneous condensates. To be general, the sonic 
horizons we consider can be due either to the velocity 
flow w(a;), or to a varying sound speed c(a;), or to a com- 
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bination of these two. We shall later see that some fluxes 
properties (such as the backscattering mixing left- and 
right-moving phonons) crucially depend on the particu- 
lar combination which is realized. 

In Sec. mil we analyze the mode equation, we identify 
the combinations of solutions which describe initial and 
final asymptotic phonons, and relate them by the afore- 
mentioned 3x3 Bogoliubov transformation. In the next 
section, we show how it governs both local and nonlo- 
cal observables. This is carried out twice, without and 
with an initial temperature. In appendices, we consider 
the scattering of coherent states, which links t he p revi- 
ous analysis with hydrodynamical experiments [2l|, and 
white holes. 

In Sec. |Vl we numerically solve the mode equation and 
obtain the spectral properties of the emitted phonons. 
We consider both weakly and strongly dispersive regimes. 
The transition from one to the other characterizes the 
validity range of the gravitational analogy. 



II. THE SETTINGS 
A. Dilute gases 

We give here the basic ingredients which describe the 
condensates and their linearized perturbations. The 
reader is referred to the review for more details. 

In a second quantized formalism, the set of atoms is 
described by a field operator '^{t, x) which annihilates an 
atom at t, x, and which obeys the equal time commutator 



[vl/(t,x),vl/t(t,x')]=53(x-x'). 



(1) 



The time evolution of ^' is given by the Heisenberg equa- 
tion 



(2) 



where the Hamiltonian is 



H = I d?x 



Am 



(3) 

In this expression, m is the mass of the atoms, V the 
external potential, and g the effective coupling constant 
which describes the scattering of atoms by a local term. 

At sufficiently low temperature, of the order of 300 nK 
for 10'' atoms, a large fraction of the atoms condense in a 
common state. To separate this state from its perturba- 
tions, the field operator is decomposed into a c-number 
wave describing the condensed atoms, ^'o, and a fiuctu- 
ation: 



m = m 







(4) 



In the mean field approximation, satisfies the Gross- 
Pitaevskii equation 



ihdt-^o = [T + V + gpo] *o, 



(5) 



where the kinetic operator is T = — fi,^V^/2m. This 
equation guarantees that po = l^oP obeys the continuity 
equation: 9tpo + div(/9ov) = 0, where v is the condensate 
velocity. 



B. Stationary condensates 

In the general case, V and g depend on both x and 
t. In the body of this work we shall only consider sta- 
tionary cases. In Appendix [XI the time-dependent case is 
presented. Before proceeding, let us make clear that sta- 
tionarity means here that there is a Galilean frame (that 
needs not coincide with the lab frame) in which V, g and 
therefore po only depend on x. From now on we work in 
this "preferred" frame. 

In this frame, the condensate wave function has the 
form 



*o(t,x) 



-ifit/h 



Vpo(x) ( 



,iWo(x) 



(6) 



where po(x) gives the (mean) density of condensed atoms, 
/i is the chemical potential, and ko(a;) — d^Wo is the 
condensate wave vector. 

In the sequel, we shall work with one dimensional con- 
densates. This means that the transverse excitations 
have energies much higher than the longitudinal ones, 
and than the interaction energy gpQ. The transverse ex- 
citations can therefore be ignored at sufficiently low ener- 
gies. This simplifying hypothesis can be relaxed without 
significantly modifying the forthcoming treatment. For 
one dimensional stationary condensates, the continuity 
equation reduces to 



po{x)v{x) — est. 



(7) 



where v{x) — hko{x)/m is the velocity of the condensate, 
and where x is the longitudinal coordinate. Plugging 
Eq. (O in Eq. (O gives 



^ +Po ^Po +^(a;) +.g(a;)po(a;) 



(8) 



Because of Eq. (O, the nonhomogeneity of the back- 
ground can be characterized by only two functions. We 
shall use v(x) and the x-dependent speed of sound 



(?{x) = g(x) 



Po(a;) 



(9) 



because the equation for the relative fluctuations only 
depend on these functions. 



C. Bogoliubov-de Gennes equation 

We show that the relative fluctuations obey, at the lin- 
ear order, a fourth order equation which does not involve 
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the external potential. Inserting Eq. ([4]) in Eq. ([2]), and 

linearizing the equation in 0, one gets 



ihdt^ =[T+V + 2g\^o\^] (t> + 5^-0 



2^t 



(10) 



Given the structure of this equation and Eq. ([5]) , we found 
that it is mathematically more convenient to work with 
the relative fluctuation 4> defined by 

* = *o (! + </>)■ (11) 
Then using Eqs. (|6l9p . one gets 

ih{dt+vd^)(l)^Tp(l) + mc^[(j) + (l)'^], (12) 
where we have introduced the "dressed" kinetic operator 



2rnpo 



V 1 

dxPo dx = -ti^ — dx-dx, (13) 
2m V 



which takes into account the nonhomogcncous character 
of the condensate density. The last expression is obtained 
using Eq. ([7]) , and is valid for stationary condensates only. 

In this last case, the field operator can be written as a 
superposition of the form 

0^ (i, x) = e-^"* 0^ (x) + al (e~»"* {x)) *, (14) 

where Uui , aj^ are phonon annihilation and creation op- 
erators. Inserting Eq. (fT4|) in Eq. (fT2|) . and taking the 
commutator with both a^j and a^, yields 

[h (uj + iv dx) -Tp- mc^] (j)^ = mc^ ip^, ^^^^ 
\—h{uj-\-ivdx) — Tp — mc^] Lp^^ = mc? 

It is instructive to obtain a single equation for (p^^ by 
eliminating (f^. This can be done by dividing the first 
line by and acting on the resulting equation with the 
operator between brackets in the second line. After some 
manipulation, we obtain 



[h [uj + iv dx) + Tp] \ [-h [lo + iv dx) + Tp] 



-h^vdx-dx I = 0. 
V J 



(16) 



This equation (or equivalently Eq. (jlSp ) is valid for all 
stationary ID condensates. It contains no approxima- 
tion besides those involved in the BdG equation [l^ [l^ . 
In Appendix [B] we analyze its properties as well as its 
relations with other dispersive models. 



D. Near horizon trajectories and background 
profiles 

A sonic horizon can be obtained in two cases depending 
on whether it is c-\-v or c—v that crosses zero. Assuming 
it is c -I- V, the condensate flows to the left, i.e., v < Q. 



We further require that c + v smoothly crosses zero, i.e., 
the following expansion is valid in a finite range: 



c + v = nx + 0{x^). 



(17) 



We set to a; = the location of the sonic horizon, and 
we call the near horizon region^ the range of x where the 
neglect of nonlinear terms in Eq. p7p is valid. 

To verify that this profile gives rise to a black hole hori- 
zon, we can either refer to the gravitational analogy 0], 
or directly analyze the characteristics of Eq. (fT6|) . Adopt- 
ing the second approach, we perform a long wavelength 
approximation to Eq. p6p . i.e., we drop the Tp terms, 
and consider the WKB solutions of the resulting equa- 
tion. In this approximation, right (u) and left (w) mov- 
ing solutions (with respect to the condensate) decouple 
and are governed by cc-dependent momenta fc", and fc^ 
respectively. These are solutions to 



UJ — k^v{x) — ±fc(j c(a;), 



(18) 



where the -I- (— ) sign governs fc^ (fc^), and where c{x) > 
0. To get the space-time trajectories we consider the 
Hamilton- Jacobi equations {dx/dt = duj/dk, dk/dt — 
—duj/dx). These trajectories give the locus of construc- 
tive interference when considering wave-packets. In the 
near horizon region, for the right movers, Eq. (jlSp gives 



k^ko e"''* 
X = xo e"*. 



(19) 



By definition, whenever such equations are obtained with 
K > 0, one is dealing with a black hole horizon. When 
instead k < 0, the structure of the trajectories is that of 
a white hole. The two cases are related by a time reversal 
symmetry, see Appendix iDl 

In both cases the relevant quantity is the "decay rate" 
K, given by the gradient 



d{c + v) 



dx 



(20) 



horizon 



evaluated at the sonic horizon c + v = 0. (In the gen- 
eral relativistic jargon it is called the "surface gravity" 
(when multiplied by c). It plays a crucial role in the laws 
of black hole thermodynamics [13, [2^, and it fixes the 
temperature of black hole radiation when using relativis- 
tic fields [l^l-) One should also point out that the left 
movers, the ii-modes, are hardly affected by the horizon 
since —c + v^2v — Kxis approximately constant in the 
near horizon region. 

In the sequel, we shall work with 



if/ \kx\ \ 

c{x) + v[x) — cqD X sign(a;) tanh" 

[\DcaJ 



(21) 



The parameter 1 > I? > determines the size of the 
near horizon region, namely |Ax| — Dcq/k. As we shall 
see, it plays a critical role in the deviations with re- 
spect to the standard relativistic fluxes. The power n 
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Figure 1: Upper plot: shape of the profile (c + v)/co as a 
function of kx/cq, for D — 0.5 (solid line) and D = 1 (dashed 
line). Lower plot: separate profiles c{x)/co and v{x)/co, for 
D — 0.5 and n — 2. The solid lines correspond to g = 0.3 and 
the dashed lines to q — 0.7. Both pairs of profiles give rise to 
the same function c + v, represented by the solid line in the 
upper plot. 

controls the sharpness of the transition from the near 
horizon behavior to the asymptotic flat regions on either 
side. For n —^ oo, the transition becomes sharp. As 
discussed in , sharp transitions give rise to nona- 

diabatic effects which produce superimposed oscillations 
on the fluxes. In t his p aper, we shall work with n equal to 
2, and we refer to [17| for more details about this aspect. 

Given that the left-moving u-modes "see" the combi- 
nation c — V, and are coupled to the u-modes, we need 
to fix both V and c, and not only the combination of 
Eq. (j2ip . To this end, we introduce a new parameter q 
which specifics how c + v is shared between c and v. 

c{x) = Co + (1 - g) [c{x) + v{x)] , ^^^^ 
v{x) = -Co + q [c{x) + v{x)] . 

For (7=1 (resp. q — 0) the_hole is purely due to the 
gradient oi v (resp. c). In the analysis was carried 
out with 9 = 0, whereas in the recent experiment j25j . one 
finds g ~ 0.7. In our numerical analysis we shall see that 
the value of q has an important impact on the spectrum 
of the Hawking radiation. The influence of D and q on 
the functions c{x) and v{x) is illustrated in Fig.[TJ 

III. THEORETICAL ANALYSIS 

To prepare the numerical analysis of Eq. (|15p . it is 
worth studying the modes {(j>uj, ^u), and the Bogoliubov 
transformation relating the asymptotic in and out mode 
bases. 

We remind the reader that for a relativistic 2D massless 
fleld, the Bogoliubov transformation induced by propa- 



Figure 2: Graphical resolution of the dispersion relation 
Eq. ((23ll, for g = 1 (left plot), and g = (rightplot). The 
straight lines represent — v±k and the curves represent 
Q.±{k). The solutions k(uj) are given by the abscissa of their 
intersections. 

gating the field in Eq. ((2T|) is particularly simple because 
left and right moving modes decouple, see e.g. When 
dealing with Eq. (|16p . one faces two novelties. First, at 
fixed uj, four independent solutions now exist, and sec- 
ondly, the left-right decoupling is no longer exact (even in 
the dispersionless limit m oo) . Similar features have 
been already confronted in (sl. [iTj. However, the present 
case is more general, since both v and c vary. It is also 
more complicated, because 4> and Eq. (jA6|) are complex, 
whereas the field and the mode equation in the former 
works were both real. 



A. Asymptotic solutions 

In the above profiles, c and v are asymptotically con- 
stant in the regions \kx\ 3> Dcq. In these regions, the 
general solution of Eq. (|16p is a superposition of plane 
waves e^^^ with constant amplitudes. The asymptotic 
roots k{bj) are solutions of (sec Eq. (jB3[) ) 

[u-v^kf ^cie + ^^nl{k). (23) 

For w > 0, in the subsonic region cj^ < there is one 
real root fc^ > which corresponds to a right mover and 
another real root fc^ < which describes a left mover, see 
Fig. [21 where the two extreme cases q = 1 (left plot) and 
(7 = (right plot) are represented. There is also a pair 
of complex conjugated roots (since the equation is real). 
The root with a negative (resp. positive) imaginary part 
corresponds to a growing (resp. decaying) mode to the 
right. 

In the supersonic region c_ > for lo larger than 

a critical frequency Wmax (we compute its value below) 
there are only two real roots, as in a subsonic flow. In- 
stead for Q < uj < cjmax four real roots exist. This dou- 
bling of the number of real roots is illustrated in Fig. [S] 
Thus the pair of growing and decaying modes which ex- 
isted in the subsonic flow is replaced by a couple of os- 
cillatory modes. Such a replacement is a generic feature 
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Figure 3: Graphical resolution of the dispersion relation in 
the supersonic region for three values oi uj: Q < uj < i^max, 

1^ = t^max and UJ > l^max- 

of QFT in external fields, see [l^]. In BEC it will occur 
at all horizons, for both black and white holes. 

B. Maximal frequency tJmax 
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Figure 4: Contours of constant li^max/K in the {D,X) plane, 
for g = (solid lines) and q — 1 (dashed lines). 




The maximal frequency Wmax is the value of lo where 
the two extra real roots merge. It is thus reached when 
the straight line uj — v-k is tangent to — r2_(fc), the neg- 
ative root of Eq. ([23|) evaluated in the supersonic region. 
The corresponding value of fc,nax is 



-8c2 



(24) 



Then Wmax is obtained by replacing k by fcmax in Eq. ([23 
Using Eq. and Eq. ([221), it is thus of the form 



Axf{D,q), 



(25) 



where the "healing" frequency A is related to the healing 
length computed with cq 



Co 



\/2toco ' 



(26) 



by A = ^/2cq/£,q. (The prefactor V2 has been added so 
that the quartic term in the dispersion relation Eq. (j23[) 
be equal to k'^cfj/A^.) In what follows, A (or the adimen- 
sional A = A/k) is referred to as the dispersion scale and 
is used to characterize the importance of dispersion. 

The contours of constant Wmax/^ in the (Z?, A)-plane 
are shown in Fig. 21 for g = and q = 1. One sees that 
q has only little influence on the value of Wmax/^ when 
this latter is small, but the effect becomes significant for 
larger values. The physical consequences of this shall be 
seen in Sees. IV C 21 Notice that for D ^ 1, one has 



f{D,q) cx 
smaller than A. 



This means that Wmax can be much 



C. Mode orthonormality and mode completeness 

To proceed to the canonical quantization of cf), one 
needs a mode basis which is orthonormal and complete. 
The orthonormality is defined with respect to the con- 
served scalar product on the space of the solutions of 
Eq. (jlSp. In terms of the doublet Wi = {(f>i, (pi) the scalar 
product takes the form [28j 



iWi\W2) 



dxpo{x) [(^*(^2 - </5^<y32] 



(27) 



The presence of po in this product follows from the use of 
the rescaled fluctuations defined in Eq. pT|) . Similarly, 
the equal time commutator reads 



1 



S{x-x'). 



(28) 



This a;-dependent measure induces no difficulty, and 
moreover one can get rid of it by using the non-Cartesian 
coordinate y defined by dy = dxpo{x). 

When the condensate is homogeneous the quantization 
is rather straightforward and well known [l^. However, 
when the flow c + v possesses a horizon, the situation is 
more subtle. Therefore, before considering the inhomo- 
geneous backgrounds of Eq. ((2T|) . let us quantize </> when 
c and V are constant both in space and time. The com- 
plete description we shall obtain transposes to a BEC 
that of Ref. [H. 



6 



1. Homogeneous condensates, k-representation 

In homogeneous condensates, it is appropriate to ex- 
press the field in terms of exponentials e*'^^ and cre- 
ation/destruction operators labeled by the (real) wave- 
vector k, 



Therefore, when using lo instead of k to label modes 
and operators, one should discard the growing and the 
decaying solutions of Eq. (|16p . Unlike the oscillatory 
plane waves they cannot be reached in a limiting pro- 
cedure, starting with square integrable functions. 



dk 



dk (j)k + d\ ifl 



(29) 



where 



(t)k{t.,x) = 



-iujf^ t+ikx 



■U/c, (pk{t,x) = 



^ — iuJkt-\-ikx 



(30) 

and where ojk is the positive solution of Eq. ((23)) . Using 
these expressions, and introducing Wk — {fl.,(t)l), the 
"bar" doublet associated with Wk — {4>k, fk), one verifies 
that the orthonormality conditions 



{Wk\Wk.) = -{Wk\Wk.) = 6{k-k'), 
iWk\Wk') = 0, 



(31) 



are satisfied when the amplitudes Ufc, Ofc take their stan- 
dard value [1^ , irrespectively of the (sub- or supersonic) 
value of the condensate velocity v. Explicitly, using 
Eq. nil) and Eq. jBSl, one obtains 



Dk 



DkUk 
1 



Ufc 



1, 



h\ fc2 + 



h?k^ 

Aw? 



rac 



2m 



(32) 



The fact that Dk is independent of the condensate ve- 
locity V follows from Galilean invariance, the condensate 
being homogeneous. Indeed, in the coordinates tc, Xc co- 
moving with the fiuid, related to the coordinates t, x by 
tc — Xc — X ~ vt, the 2D wave vector of compo- 
nents (wfc, k) in the t, x system has comoving components 
(rifc — ujk — vk, k) where all reference to v drops out when 
using k to label modes, as can be seen from Eq. (^5]) . 

The operators afe,a|, obey the usual bosonic commu- 
tation relation [afe,a^,] ~ 6{k — k'). The relationships 
between mode doublets Wk and these operators follow 
from 



dk = {Wk\W), 
al = -{Wk\W), 



(33) 



where the doublet operator is = {(j), 0^'). 

In this fc-representation, one can also easily verify 
that the mode basis is complete. That is, starting with 
the commutators [afc,a^,] = 5{k — k') and the doublets 
VFfc, making use of the completeness (in the sense of 
Fourier analysis) of the exponentials e*'^^ with k real from 
[— cx), oo], and using that Dk = D^k (which is a necessary 
condition) one verifies that Eq. (pS)) is satisfied, again ir- 
respectively of the value of the condensate velocity v. 



2. Nonhomogeneous condensates, ui -representation 

When considering nonhomogeneous but stationary 
fiows, one should use the conserved frequency uj to la- 
bel modes. Then one notices that the change of variable 
k LO in Eq. (j29p proceeds very differently according 
to the sub- or supersonic character of the flow. As a 
preliminary step, we consider separately sub- and super- 
sonic homogeneous fiows. A superscript u (v) shall be 
added to characterize right (left) movers with respect to 
the condensate. 

In a homogeneous subsonic flow, only two real roots of 
Eq. (USD exist: fc"(w) > and k^iuj) < 0. Given that 
dk/duj never crosses zero, one can re-express Eq. (1291) as 



0(i,x)= / duo e-^"*0^(x) + e+*'^*^i(x) , (34) 



where 



^^{x)=dicpi{x) + dl4>l{x), 

Vu>{x) = dZ (fiZix) + dl if>l{x). 



(35) 



The rescaled modes and operators are = 4>k \/ dk/du), 
and = dk\/ dk j duo . One easily verifies that the fac- 
tors ^ dkjdLo guarantee that all 6(k~ k') obtained in the 
former subsection are consistently replaced by S{lo — uo'). 
That is, the operators d^, dj^ obey [d^^, d^,] = S{lo — co'). 
Similarly the modes 0^ and 0^ are orthogonal to each 
other, and possess unit positive norm (in the sense of a 
Dirac distribution 6{lo — lo')). 

We now consider a homogeneous supersonic fiow. 
Starting again from Eq. (P^ one decomposes, as in 
Eq. (p4)) . the field as an integral over a; of a sum of 
right- and left-moving modes. When considering the left- 
moving sector in left-moving fiows, dk/diu does not cross 
zero. Therefore, as in subsonic fiows, all left-moving (pos- 
itive norm) modes can still be monotonically labeled by 
LO belonging to [0, oo]. The same is no longer true for the 
right-moving sector. In fact, the integral dk splits 
into an integral over co belonging to [0, oo] plus another 
piece over negative frequencies belonging to [— WmaxjO]. 
This last interval terminates at — Wmax where the two new 
real roots A:"(w) > merge. 

Thus, for to > Wmax, read as in Eq. ([35]) since 

one has only one (positive norm) w-root. Instead, when 
< LO < Wmax, three real u-roots exist: the continuation 
(in lo) of this positive norm one, plus two new roots with 
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negative 11, see Fig. |3l In this case, one has 
dt 



2ti 



■ (j){t, x) 



(36) 



(37) 



1 = 1,2 



In the first equation, a complex conjugate and a subscript 
—bj have been used to characterize the two new modes. 
This means that the two doublets W^t^^i = (0"^^,;, ^"^^1,1) 
have a positive norm and obey Eq. (jl5[) with a frequency 
idt = —CO < 0. It should also be noticed that the 
above operators (j)^ (x) , ip^ (x) contain both annihilation 
and creation sectors. This allows one to write (l>{t, x) as 
in Eq. ([34)) . i.e. as an integral over cu S [0, 00]. 

In metrics which contain a transition from a subsonic 
to a supersonic flow, because of the scattering on v{x), 
the modes that are bounded on one side of the hori- 
zon are generally not bounded on the other side. How- 
ever, given the linearity of Eq. (fT5|). one can always 
construct bounded modes as linear combinations of the 
above modes, requiring that the coefficient of the growing 
mode be zero [l7|. 

In fact, when lo < Wmax: three independent bounded 
combinations can be constructed and (jjuj^ instead of 
Eq. ([55)1 . now reads: 

k (x) = aZ cj^Z (x) + al ix) + all [^1^ (^)] % (38) 

and similarly for ipi^{x). In this expression, 0^, 0"^ 
and ipl^ stand for yet unspecified, normalized bounded 
modes. The particular cases of the in and out bases shall 
be defined in the next section. For lo > Wmax instead, 
there is one growing mode on each side of the horizon. 
Hence only two independent bounded linear combina- 
tions can be constructed, so that (j)^ is still decomposed 
as in Eq. ([55)1 . 



D. Bogoliubov transformation 

1. Vacuum instability and spontaneous pair production 

Since we are dealing with a stationary situation, the 
energy of the state of is constant. Hence, one might 
have thought that no phonons could possibly be sponta- 
neously emitted. This is not the case, because the vac- 
uum is unstable against the production of phonon pairs 
when two conditions are met. First, negative energy ex- 
citations must exist. We saw in the last section that it is 
the case when the flow is supersonic, for u < Wmax- 



However, this is not enough, as can be understood by 
considering a homogeneous condensate propagating in a 
frictionless translation faster than sound in the lab frame. 
One also needs spatial gradients to define unambiguously 
the unique "preferred" stationary frame, and to couple 
the so-defined negative energy excitations to positive en- 
ergy ones. When these conditions are met, the Hamilto- 
nian of has the following structure 



H = 



du;ha;{aZ^aZ + aZ^aZ). (39) 



From this expression three conclusions can be drawn. 
First, phonons with u > Wmax cannot participate in the 
pair production since no partner with the corresponding 
negative frequency exists. Second, both u and v positive 
frequency phonons participate in the vacuum instabil- 
ity. Which of the two channels contributes most depends 
on the intensity of the coupling with the negative fre- 
quency excitations. Third, the diagonalization of H in 
the sector < a; < uj^^ax is ambiguous because by a uni- 
tary (w-diagonal) Bogoliubov transformation, the above 
quadratic form is left unchanged. 

Therefore, one needs an additional physical criterion to 
remove this ambiguity. In the general (time-dependent) 
case no such criterion exists, and the notion of phonon is 
inherently ambiguous [H, |33|- However, in the present 
case of stationary asymptotically homogeneous conden- 
sates, there is no ambiguity to define two complete sets of 
modes. The first set of modes (the "in" modes) charac- 
terizes the initial phonons propagating towards the sonic 
horizon. The second set, the "out" modes, character- 
izes the asymptotic particle content of the scattered field 
configurations. Having the two sets, we can compute 
how they are related and how this relationship governs 
the decay of the vacuum. 



2. Space-time structure of in and out wave-packets 

The procedure to identify the in and out modes is stan- 
dard [2y| . One should construct "broad wave packets" , 
i.e., superpositions of stationary modes, so as to extract 
the asymptotic temporal behavior by looking at the sta- 
tionary phase condition duiS = 0. This equation is equiv- 
alent to Hamilton's equation since the mode phase S co- 
incides with the Hamilton- Jacobi action in the WKB ap- 
proximation. From this asymptotic behavior, one identi- 
fies the modes (in fact the doublets W^^ a — {4'ui,a, fuj.a)) 
that are associated with each initial (and final) one- 
phonon state. It should also be pointed out that these 
modes acquire a precise physical meaning when consid- 
ering coherent states, see App. [C] 

It is useful to visualize these modes. Let us first de- 
scribe the in mode 0"'*", i.e., the particular solution of 
Eq. (|16p which contains in the past only a right-moving 
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Figure 5: Schematic space-time representation of a wave- 
packet made out of modes (pZ'^", for u > tJmax- At early times, 
the packet is purely right-moving. It is then scattered into a 
transmitted right-moving part and a reflected left-moving one. 
The dashed line represents the initially unexcited "ancestor" 
of the late-time left-moving packet. 




1 



Figure 6: Same as Fig.[5]for uj < ojmax, when pair production 
occurs. 



3. Bogoliubov transformation 

From Eq. ([M)) one sees that for uj > tJmax one has 
only two modes with positive frequency/energy. Because 
the condensate is not translation invariant, there is some 
scattering. Therefore, in and out modes are related by 

A^uAn rp iu,out I D IV, out 

IV.in D* JAl.OUt I rr\-^ IV ,OUt ^ ^ 

The other modes, the ^^Z"'™/""*^ g^^g related among 
themselves by the same relations. Since there is no mix- 
ing of with iy9* , the conservation of the norm trivially 
implies |T;jp -I- \Ruj\'^ — 1. One is thus dealing with an 
elastic scattering between u and v modes, and T^^, are 
the transmitted and reflected amplitudes, respectively. 
This is a "trivial" transformation in the sense that there 
is no spontaneous pair production. The vacuum of these 
high frequency modes is thus stable. 

It should be pointed out that this elastic scattering is 
also found, for all modes, when the flow remains every- 
where subsonic. In fact, because of dispersion, the high 
frequency modes with oj > cjmax do not see/experience 
the presence of the sonic horizon. This can be verified 
by computing the trajectories followed by wave-packets 
centered around some lo > Wmax- 

The situation is radically different for lo < ojniax- In 
this case the mode mixing is nontrivial, and three equa- 
tions are needed to characterize the transformation 

iu,in :U,out I iQ /nout\* I r> I , „v ,out\* 

The coefficients are given by the overlap of the corre- 
sponding (normalized) in and out doublets, e.g. 

The normalization of the coefficients then immediately 
follows, e.g. the first equation (together with the corre- 
sponding one for </?"'*") gives 

+ 1^1'- 1/3-^1' = 1. (43) 



packet. We shall do it twice, both for lo > Wmax, where 
there is only some elastic scattering (see Fig.[S]), and for 
UJ < oJniax, in the presence of pair creation (see Fig. [6]). 
In both cases, initially, one only has the incoming branch 
with unit norm. At late time, when uj > Wmax, one has a 
transmitted M-mode with amplitude T^^, and a reflected 
?;-mode with amplitude R^j- For uj < Wmax, in addi- 
tion to these modes there is the negative frequency mode 
((p"'°"*)*. The description of the other in and out modes 
is obtained without difhculty. 



In this expression, the minus sign comes from negative 
norm doublets, see Eq. PT|) . 

It should be noticed that the above enlarged Bogoli- 
ubov transformation governs the general case. It applies 
indeed to any stationary situation when there is one type 
of negative frequency modes (here the u-modes (t>1^)- It 
should also be noticed that in situations with two sonic 
horizons, such as black hole white hole pairs 
the Bogoliubov transformation will be more complicated 
than Eq. (|iT|) . 
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IV. FLUXES, DENSITY FLUCTUATIONS, AND 
CORRELATIONS 

Most of the work dedicated to the Hawking effect con- 
centrates on the particle content, or the energy content, 
of the outgoing flux p, 3Q|. These two observables are 
related to of Eq. (|¥T|) . However, it was also noticed 
that there exist Einstein-Podolski-Rosen (EPR) corre- 
lations between the outgoing particles and their part- 
ners with negative frequency oj. Unlike the energy flux, 
these correlations are weighted by ao;/?^ and originate 
from interfering terms (i.e., nondiagonal in the occupa- 
tion number). In the case of gravitational black holes, 
little attention has been given to the late time correla- 
tions because they are hidden for the external observers 
since the partners are trapped inside the horizon [26j . 
Nevertheless these correlations have well-defined prop- 
erties [S^l- Moreover, they extend to the past and can 
be revealed by sending quanta to stimulate the emission 
process [s^. IsEI, l36j . 

With the advent of acoustic black holes, the sit- 
uation completely changes because one has access to 
both regions, and can therefore probe the EPR corre- 
lations. In this respect, acoustic black holes are simi- 
lar to what is found in the homogeneous time-dependent 
EEC [13, m, m and inflationary cosmology. In 
that case, the primordial fluctuations, seeds of the galaxy 
clusters and of the temperature anisotropies in the Cos- 
mic Microwave Background, also result from pair cre- 
ation |41|. Moreover, their correlations possess a well- 
defined space-time structure |20] which, on the one hand, 
affects today's observables because both members are 
within our Hubble patch, and on the other hand, is very 
similar to that found in time-dependent BEC [H, . 

We start our analysis with the 3 occupation numbers, 
with and without an initial temperature. The first case 
is similar to that of but many new features arise 
because of the different structure of the wave equation 
and because both c and v vary. The second case is com- 
pletely new, as are the next sections where we relate the 
coefficients of Eq. (PT|) to both local and nonlocal density 
fluctuations. 



A. Occupation numbers in the initial vacuum 

Let us first consider the in vacuum |0i„), i.e., the state 
annihilated by the destruction operators a"'™, a^'*" and 
a"'™, defined by the in modes through Eq. (|33p. For grav- 
itational black holes, this is the physically relevant state 
(for outgoing u-configurations) after a few e-folding times 
At = In fact, because of the exponential redshift 

effect associated with the near horizon propagation, see 
Eq. , the transient effects due to infalling quanta are 
exponentially rapidly washed out. Therefore, the choice 
of the initial distribution of quanta does not affect the 
stationary properties of the outgoing flux. This is no 
longer necessarily true for acoustic black holes, because 



the dispersion limits the number of e-folding times during 
which Eq. applies [H,!!!]. Thus, one should analyze 
each case to see if the vacuum state provides a reliable 
approximation. 

Assuming this is the case, the mean occupation num- 
bers are 





= (o.„ 


out \ - u. out 


|0,„) 




^vac, V 


= {0^n 


1-1), out t - u, out 


|o,„) 




"■-UJ 


= {0^n 

— n"^" 

"uj 


\ "U, out t - li. out 


|0.„) 





These expressions follow when using Eq. (|4ip and 
Eq. p3p to express the out operators in terms of in ones. 

When compared with the simpler case where u and v 
modes decouple the main novelty is that t^-quanta are 
also produced (to the left of the horizon, in the "inside" 
region when using the gravitational analogy). Their oc- 
cupation number is n^^''". Because of this, n™'', the 
numbers of u phonons emitted to the right is always 
smaller than that emitted to the left: n™^. In fact, 
n™^ = n^" + n^"' " tells us that in addition to the usual 
"Hawking" channel, there is a new channel where one of 
the partner is a w-quantum. When the latter is negligible, 
for \Bui\'^ <C one recovers the simpler case where 



B. Occupation numbers from an initial thermal 
state 

In a BEC, in realistic situations, there will always be 
some residual temperature. The order of magnitude of 
this temperature is given by the chemical potential fi of 
Eq. ([8]), and thus inversely proportional to the healing 
length ^. Moreover, the characteristic wavelength of the 
condensate inhomogeneity ~ c/k will generally be larger 
than ^. Therefore, since Hawking temperature is ksTH — 
hK/2Tr, the initial distribution of phonons could well hide 
the Hawking effect. 

For concreteness, to characterize the initial state, 
we assume that the three distributions have the same 
temperature. Using fi*", the initial value of the co- 
moving frequency, the three initial occupation numbers 
n™,n™'",n™„ are 

-vn,a ^ |^g/3T n'-'-iuj) -iy\ (45) 

where Pt is related to the (initial) temperature by 
ksTin — h/ (3t- This choice means that, before the scat- 
tering in the near horizon region, the temperature mea- 
sured in the frame comoving with the fluid is the same 
for all modes. 

However, because of the scattering, modes sharing the 
same constant frequency w mix. Hence we must use uj 
to characterize the initial distributions. For each type 
of modes, we thus need to express 57*" in terms of to. 
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This explains the presence of the index a in the comoving 
frequency in the above equation. Using Eq. the 
values corresponding to the three in modes of Eq. ((4T|l 
are 

n'^'-^iio) =uj-v+k''{Lo), (46) 

where the three roots fc"(w) > 0, fc"(w) < 0, k"{-uj) > 0, 
are clearly seen in Fig. O 

Given the initial occupation numbers, Eq. (|4T|) fixes 
the final ones to be 

-fin,. ^ ^»,- + |ij2 

+ |B^P(l + n-+ (47) 
+ |i3^|2(l + n!_"^+n»^-). 



distinguish left- from right-moving phonons, it could be 
possible to detect Hawking radiation, even when the ini- 
tial temperature is larger than Hawking temperature. In 
Sec. IV Dl we study numerically realistic situations and 
confirm this possibility. In the case one cannot distin- 
guish left- from right- moving phonons, the situation is 
much worse. The dominant noise term would be n™''', 
and an initial temperature larger than Hawking temper- 
ature would hide Hawking radiation. 

Before proceeding to the numerical analysis, we study 
density fluctuations and nonlocal density correlations, 
firstly because phonon occupation numbers are not di- 
rectly measurable (see however [l^]) whereas density fluc- 
tuations are, and secondly because nonlocal correlations 
are amplified by initial thermal distributions instead of 
being smeared by them. 

The reader interested in the spectral properties can 
read Sec. |V] first and go back afterwards to the next sec- 
tion. 



D. Density fluctuations 



The interpretation of these equations is clear. The first 
term is the corresponding initial occupation number. 
Then for the first two equations, the second term is due 
to the elastic scattering between u and v modes which 
adds (or subtract) particles according to the strength of 
the scattering, whereas the last term is due to the in- 
duced emission which involves both the partner's initial 
occupation number n™^ and that of the species itself. 
In the third line instead, one has two induced emission 
terms because there are two "pair creation" channels and 
no "elastic" channel. 



C. How to get rid of thermal noises? 

From the first equation in (j47|) , one can easily imagine 
that the Hawking radiation (HR), i.e., the spontaneous 
creation of pairs weighted by l/J^jP given in the first line 
of Eq. ([33]), might be hidden by the presence of initial 
distributions. See also [l^ for a discussion of the conse- 
quences of three-phonon interactions. 

There is yet another difficulty which can complicate 
the detection of HR, namely the possibility to distin- 
guish right- from left-moving phonons. In the case one 
can, detecting HR requires that be larger than, or 

at least of the same order as, both n™ and 
The first condition could be satisfied because the initial 
distribution of u modes can be significantly redshifted. 
By this we mean that one can have i7™'"(±ijj) ^ cu and 
thus possibly r2™'"(±a;)/Ti„ ^ uj/Th, although a priori 
Tin ^ Th- The second condition might be more prob- 
lematic because the v modes are hardly redshifted. Nev- 
ertheless, the inequality lA^^pn^'" < \(3uj\'^ could also 
be satisfied because, as we shall see, in certain cases the 
u — v mixing is very small. Therefore, when one is able to 



Given that </> is a complex field, there are several ways 
to characterize the fluctuations in a BEG: either through 
the density correlation function which is governed by 
Re (j), see Eq. (jB5|) , or through the phase correlations 
governed by Im0, or even through the crossed phase- 
density correlations. In what follows we only discuss the 
density-density correlations as the extension to the two 
other types is easily made. 

To simplify the forthcoming expressions we introduce 
the field operator % = ^ -I- In stationary cases, it can 
be decomposed as in Eq. ([M]) . and in terms of the same 
operators as those of Eq. ([55]) . The only change is that 
the wave functions 0°, (p^j are all replaced by 

xZ{^) = rj^)+^z{^)- (48) 

Since this correspondence applies to both the in and out 
sets, and since both and (p^ obey transformation of 
Eq. ([311), the three initial x™'" are also related to the 
three final x™*'" by Eq. ([4T|) . Therefore, even though 
X is not a canonical field, as it does not obey canonical 
commutators, it can be treated as a genuine quantum 
field when computing correlation functions. 

The statistical properties of the density fluctuations 
encoded in a given state are characterized by the anti- 
commutator 

/OO 
duje-"^^*^''^ G'^\x,x'), (49) 
-OO 

where p'" is the initial density matrix, since we work 
in the Heisenberg representation. In the second line we 
passed to a Fourier transform since we assumed that both 
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the condensate and the state are stationary (in the "pre- 
ferred" frame). 

When the initial state p™ is incoherent, using the in 
basis to express Xi only three terms having the same 
structure are obtained 

GL"(x,x') = (< + 1/2) xL"'"(^) [xL"'"(^')]* 

+ {fi}-'- + 1/2) xr"(^) m 

+ (nL-^ + l/2) [x"r(^)]*x!":''(^')- 

This expression is valid for a; > 0; for w < 0, one has 
G'^^{x,x') = [G'^{x,x')]* since G'"(<, x; t', x') is real. In 
the above equation, the initial occupation numbers are 
given by 



X 5'^ = Tr[p' 



(51) 



Due to the incoherence of p'", only the diagonal terms 
remain. Additional terms would be obtained if the state 
p™ contained correlations among the initial configura- 
tions. In what follows we assume it does not (see however 
App.p. 

Because of the scattering near the sonic horizon, 
G™{x,x') has a rather complicated structure, as can be 
seen by decomposing the in modes into out ones. In 
fact G™{x,x') encodes both local observables related to 
the occupation numbers of Eq. (I47|) , and non-local ones 
governed by correlators such as Tr[p'" a°"*''a°"*'-'] with 
i ^ i- 



E. Coincidence point limit 

Far from the horizon so that the scattering is com- 
pleted, i.e., when |a:| ^ Dcq/ one should express the in 
modes in terms of their asymptotic plane wave content 
~ e^^'^^ . To ease the reading, we shall call these asymp- 
totic contributions by the corresponding wave, x™' ° or 
^out,a^ for which the amplitude of this contribution is 
unity, and, to avoid any misinterpretation, we shall add 
an upper index "as" to make clear that only the unit, 
plane wave contribution should be kept. In terms of 
these, using Eq. (|47p . in the asymptotic right region, one 
gets 

GL"(a:,a:)^(n5"'« + l/2)x |xr'"'"1' 
+ (n^'^ + l/2)x Ixr^'^l', 

whereas, on the left, the "power" is asymptotically equal 
to 

G^^{x,x) ^ (n^^'^' + 1/2) X \xZ"'-^'^r 

+ (<■" + 1/2) X Ix^'^'^^f 
+ (n!ir + l/2)x|x!":"'^T- 



(53) 



Since the point x = x' lives in an asymptotic re- 
gion where the condensate is homogeneous, the norm 



of the asymptotic modes xSf is x independent. No- 
tice also that we discarded all oscillatory terms, such as 
^out,u,iis^^in,v,asy ^ ^ because they rapidly 

oscillate as x — > cxd, and thus drop out when averaging 
over Lo. 

The interpretation of the above equations is clear. 
When working in the state p'", the two asymptotic val- 
ues of the correlator are the sum of the contributions of 
the waves that have been scattered, governed by the final 
occupation numbers, 



fin, i 



Tr[p' 



I ^ out^ i f ^ out, 2] 



(54) 



and of the waves which have not propagated through 
the horizon region, and whose occupation numbers are 
the initial ones given in Eq. ((5T|) . From Eq. (|52p one 
clearly sees that the second term (the v contribution) 
will hide the Hawking radiation whenever fi'"' " is much 
larger than n^"' which is the case in realistic situations, 
as we shall see in Sec. IVD 31 

It is therefore also of interest to compute 



T{t,x-t\x') = —(9,, -9,)Tr [r^Ht',x')^{t,x)\ 
2po 

dwe-"^^'-'"^ T^{x,x'). (55) 



In the coincidence point limit {t^x) ~ [t' ,x'), hp^T /m is 
the atom flux at {t^x). It is the sum of p^v = p^hk^/m, 
the flux of the condensed atoms, and of the integral over 
uj of !Fi^{x,x). In the asymptotic right region, using 
Eq. HID, for > 0, one has = T^' + J^™" where 



= fco X [n^ \(pj ' {x)\ +n 

^com ;^tin j^u, out \±u, out, as f 2 



in,v 



I iv.tn.as / \ I z 



(x)|^ + <"'''fc^'™IC"'"^WI^ 

(56) 



For UJ < 0, one has the same expressions with opposite 
signs, with (p^^ replaced by (pui, and n^^ replaced by n^^-l-l. 
T'^'^ arises from the uncondensed atoms and is due to the 

UJ 

dragging of the background, while J^^°™ is the atom flux 
measured in the frame comoving with the condensate. As 
in the case of the density fluctuations, the term arising 
from the initial distribution of w-phonons largely domi- 
nates. However, since fc^'™ < 0, is related to the 
difference of the terms appearing in G]^{x,x), up to k- 
dependent factors (different for each term). Thus, if one 
has access to both the atom flux and the density fluctu- 
ations in the right asymptotic region, there is a greater 
hope to have access to n^". 

Note that any other local observable constructed out of 
two fields, like the depletion, will have the same structure 
in the right asymptotic region and will thus suffer from 
the same limitations, namely it will be largely dominated 
by the contribution of n™'". This reinforces the interest 
to consider nonlocal observables. 
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F. Long distance correlations. 

1. Late time entanglement 

We now study the long distance correlations, for \x — 
x'\ ^ Dco/k > ^0- In this case G™{x, x') displays a rich 
structure. Considering a; > 0, x' < 0, one gets a priori 
eight terms, since one has two asymptotic modes on 
the right, and four on the left, see Eq. (155)) . However, 
many of them drop out because they destructively in- 
terfere upon integrating over w, see the discussion after 
Eq. l|53|) . Hence, when the initial state p™ is incoher- 
ent, no (long distance) correlations among in modes exist. 
Similarly, all terms mixing in and out modes will de- 
structively interfere. In fact, only out phonons are entan- 
gled by the scattering in the near horizon region. Hence 
only correlations among asymptotic out modes will con- 
tribute. Given that in the subsonic region, there is only 
one asymptotic out mode, -^2"*^ and in the super- 
sonic region, two such modes exist, for x > 0, a;' < 0, one 
has 

When both x and x' are taken in the subsonic region, no 
long distance correlations can develop since x^"*' ^ 
the only asymptotic out mode. On the contrary, when 
both X and x' are negative, in the supersonic region, the 
two asymptotic modes are entangled and this will show 
up in 

G'^{x,x') = Cxr''^'"^WX-"j'"'''(2;')- (58) 

Given the 3x3 character of Eq. (|4ip. three types of late 
time correlations could have been expected since three 
different couples of out modes can be formed. 

To compute A^i , and C^^ one can either express the 
r.h.s. of Eq. (|50p in terms of out modes and identify 
the coefficients multipying the corresponding couple of 
^mit,as ^ or equally start with Eq. (|49p and decompose the 
field X using the out basis. Adopting the second, more 
rapid method, we get 

-Tr[p'"a°"*'"a°"*'''t], 

6^ = Tr[p-ar'"a°'L*'"], (59) 
C = Tr[p-ar'''a°'j'1, 

where we have used the fact that the a°"* operators com- 
mute in each product. Unlike the coincidence point limit 
of G^^{x,x') which is governed by diagonal terms, see 
Eqs. ([5^1551) . the long distance correlations are governed 
by terms which are nondiagonal in occupation number. 
This is exactly as in homogeneous situations [10, and 
in fact will always be found when the parametric am- 
plification (or the scattering) conserves a quantity, the 



frequency oj here, the spatial wave-vector k in homoge- 
neous cases. A straightforward calculation gives 

= Al + ^ A^ aZ* + (n™ + 1) B^, 

(60) 

= fC A. + " a^i i?: + ("1". + 1) K a_.. 

In the vacuum, the residual correlations all involve the 
negative frequency modes because these enter in both 
pair creation channels. 

At this point, an important observation should be 
made. Because the initial distributions n*"' only appear 
in factors multiplying terms already present in the in vac- 
uum, the long distance correlations will not be erased by 
the presence of initial quanta. In fact an initial tempera- 
ture will in general amplify the long distance correlations 
induced by interactions. ^ It was numerically observed in 
[l^ . as was the fact that the A^j coefficient, arising from 
the product a"a^^, is nonzero even in a subsonic flow. 

We finally remind the reader that the above correla- 
tions were obtained using the BdG equation which 
neglects phonon interactions [T^. One might therefore 
worry that the entanglement is reduced upon taking into 
account such interactions. We refer to Refs. [i^, El, US] 
for an analysis of this point in a cosmological context. In 
brief, the weakness of the nonlinearities guarantees that 
the entanglement is hardly reduced. 



2. Spatial structure of long distance correlations 

At fixed u; no spatial structure emerges from Eqs. (|571 
[55| . To get the spatial properties of the correlations, one 
needs to take the inverse Fourier transform. Then, as it 
is the case when considering wave-packets, constructive 
interferences will develop along the characteristics of the 
mode equation. To ease the reading of the forthcoming 
expression we found convenient to return to Eq. (|49p and 
to re-introduce t and t' . 

As usual, the constructive interference condition gives 
the stationary phase condition di^S = 0. In the present 
case, using a WKB approximation for the modes xZi^) ^ 
exp(i dykZ{y)), the three phases of the terms weighted 



^ We added "in general" because the terms in the r.h.s. of Eq. (I6UI I 
have no definite sign. Hence, an increase in the temperature 
could lower, in some (rare) cases, the resulting amplitudes of the 
l.h.s. of Eq. I I60I I. It is interesting to point out that when the 
u — V mixing is small, i.e. when the A and B coefficients can be 
set to zero, only two terms in B remain (the first and the third 
one), and that these will always add up, thanks to the unitarity 
of the reduced 2x2 Bogoliubov transformation. 

^ Because of the orthogonality of the scattering matrix of Eq. I|40| l 
it would be zero unless n™' We thank C. Mayoral and 

A. Fabbri for drawing our attention to this particular point. 
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by Auj ; Bui and C^j are respectively 

SA{t, t', X, x'- u) ^ -w{t -t')+ i dykliy) 



where 



c^y^i^.fe) +arg(ln^^), 
SB{t.t',x,x'-uj)^~u;{t-t')+ ( dykZiy) 

J Z 

+ 1^ dykliy) +mg{\nBu,), 
Scit, t\ X, x'- uj) = -uj{t -t')+ f dykliy) 



(61) 



+ / dyfc"^(y) + arg(lnO, 

J z 

where z is an arbitrary location where the absolute (un- 
observable) phase of the is fixed. The stationary phase 
condition gives 



{t-t') 



dyduikZiy) 



dyduikKy) 



+duj arg(lny4. 



(62) 



for the first line, and similar equations for the second 
and third lines. It seems a priori that the choice of z 
matters. However this is not the case because z enters 
in Au) in such a way that a change of z leaves the r.h.s. 
of the equation unchanged. (This is because the coeffi- 
cients of Eq. (|59p contain the operators a\a which are 
"contravariant" with respect to a phase shift of the corre- 
sponding mode.) The physically meaningful phase that 
comes out of these expressions has the role of fixing the 
location where the interactions occur, [l^ 

This result becomes exact in the limit where x and 
x' are taken far away from the scattering zone, and it 
amounts to put z = and to to drop the last term in 
Eq. (j62p . and similarly for the equations involving B^j and 
Cuj- There could be some finite phase shift with respect 
to these WKB estimates, but these do not change with x 
and x' , and hence give subdominant effects in the large 
X limit. In this limit, the properties of the correlation 
pattern derived from Eq. (j62[) are thus independent of 
both the norm and the phase of the coefficient Au, ■ Hence 
the same long distance space-time pattern will be found 
both in the limit k — > and k oo, i.e. in regimes which 
do not give the standard Hawking radiation. Therefore 
if Hawking radiation implies the pattern, the converse is 
not true (when defining Hawking radiation as the near 
thermal radiation associated with a finite surface gravity 
«)• 

In the large distance limit, when putting t ~ t' , the 
stationary phase condition applied to Eq. (pT|) gives, for 
the Aui , Bui and Cuj terms respectively. 



At^'-^ix) 
At'''\x) 



At 
At 



HJ. u 



HJ, u 



At 



HJ, a 



ix) 



dydu^kZiy) 



(64) 



is the time it takes the a-type phonon of frequency oj to 
propagate from a; = to x, in virtue of the Hamilton- 
Jacobi equation determining the group velocity 



v'iA^) = {duikzr\ 



(65) 



For each type of correlations, we see that the locus of 
constructive interference of, say, x' given x, is given by 
the value of x' reached at the same lapse time it takes 
the partner to reach x, both phonons starting their jour- 
ney near the horizon x ^ Q. In the large x limit, the 
dominant contribution comes from the uniform motion 
outside the horizon region. Thus the above three con- 
structive interferences occur at locations x, x' related by 



(66) 



Using Eq. ((23|) . the asymptotic group velocities are, for 

X +O0, 

z;«,(l.) = d^n^ +v+^ elk + v+, (67) 

for the right movers in the subsonic region, the "Hawking 
quanta"; and, for x —oo, 



vi^iio) = -dkn- + v- = -c^_k 



1 -f fc^e 

f7_(fc) 



+ V-, 



1 + k-^e 

~~njk) 



+ V-, (68) 



(63) 



for the V quanta and the u partners in the supersonic 
region. In the dispersionless regime, for fc^± <C 1, these 
group velocities are independent of uj and respectively 
equal to C-|_ -|- in Eq. (j67p . and — c_ + f-, c_ + t;_ 
for Eq. ((68l) . Therefore, in this regime, the correla- 
tion pattern is independent of w, and will show up in 
G'"(i,x;t' = t,x') of Eq. 

These three types of correlations have been (numeri- 
cally) observed in [l3|. They have been also correctly 
interpreted save for two aspects. First, the C^j branch 
has been attributed to the "partial elastic scattering" of 
right movers of positive frequency. If this explanation ap- 
plies to the first term, it does not to the last two which 
involve pair creation B^j coefficients mixing and X-uj ■ 
Second, the Aui is claimed to "originate from thermal ef- 
fects" . From the first line of Eq. ([50]) we see that it is 
indeed amplified by thermal effects, but, because of the 
third term, it is already present in the vacuum. Being 
quadratic in frequency mixing coefficients, it is too weak 
to be easily seen in the numerical simulations. However 
using a black and white print of the arXiv version of 
p^ . these correlations are quite visible in the last two 
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plots in Figure 2, see Figure 6 for their orientation in 
the x^x' plane. Upon contacting the authors, they agree 
that these correlations are indeed present at zero temper- 
ature. When preparing the revised version of this paper, 
we became aware of 48] which agrees with this point, 
as with the rest of our analysis, and which also contains 
additional interesting results. 

Concerning the relative amplitude of A^n and C^j 
it should be noticed that, in general, there is no clear 
ordering. Instead, when the u — v mixing is weak, the 
Bogoliubov coefficients A^j, Bui are much smaller than 
f3ui (see Sec. IVC4I) and therefore, B^i is the largest. 
In this regime, one recovers the properties of the B^, 
branch relating the u modes across the horizon that have 
been known for a while. In the context of gravitational 
black holes, using relativistic fields, they can be found in 
[26,. ,33,] . As of acoustic black holes, it was understood 
in Q that the late time behavior of these correlations 
is essentially unaffected by dispersive effects,^ see also 
[43I H^l . Finally, that the Bu, correlations determine (in 
the hydrodynamical limit) the lon g di stance density cor- 
relations in BECs was stressed in [ij]. 

In App. [C]we present an alternative and simple way to 
characterize the correlation pattern. It consists in send- 
ing classical waves -described by highly excited coher- 
ent states- towards the horizon. This approach is worth 
considering because it allows to relate the above study 
of Eq. (j49[) both to the experiments in hydrodynamics 
described in t21|], and to the wave packet analysis of 



SPECTRAL PROPERTIES: 
RESULTS 



NUMERICAL 



A. Numerical procedure 



The numerical procedure we used was elaborated from 
that of [l3|. The wave equation of that reference was 
replaced by the BdG equation Eq. , where the veloc- 
ity V and the sound speed c both vary, see Eq. ((22|) . The 
extraction of the Bogoliubov coefficients from the numer- 
ical solutions to Eq. (fT5|) took into account the specific 
normalization of the modes {(puj, ^u:)- For a presentation 
of the procedure itself, we refer to |17l |. 



^ In fact the correlation pattern is as robust as the occupation 
number n^i = and the stationarity against introducing dis- 

persion. This is because z^j = P^jo-^, the complex squeezing pa- 
rameter which determines the final quantum state |2a | , is hardly 
affected. Its norm determines n^^ whereas its phase determines 
the correlation pattern. The density (in oj) of modes subject to 
the Bogoliubov transformation ensures the stationarity. 



B. Value of the parameters 

Let us determine the number of free parameters and 
their realistic ranges."* A typical value for the average 
sound speed Cq is 



Co = 0.15 cm • s 



(69) 



Assuming that the condensate is made out of ^^Rb, the 
mass of the atoms is 



1.5 X 10 



This yields the healing length. 



-25 



kg. 



Co 



\/2, 



Co 



%/2; 



mco 



A 



3.3 X 10"^ cm. 



(70) 



(71) 



The distance over which the variation of the sound 
speed and flow velocity takes place, that is, the distance 
separating the asymptotic regions where these speeds are 
constant, cannot be smaller than a few healing lengths. 
To reduce the number of free parameters, we assume in 
this work that it is of the order of lO^o (see however 
Sec. IV C"6|) . Then, given our parameterization Eq. (|2ip . 
the gradient at the horizon is 



co^ 
5^0 ■ 



(72) 



With Eq. (|7T|). this yields a relationship between D 
and A = A/k: 



(73) 



as well as an expression for the Hawking temperature: 



Th 



Hk 



D X 1.1 nK, 



(74) 



where we have used Eq. ([72|) and the numerical values 
above. We choose as free parameters D and q. Thus 
Eq. jTSl) fixes A. 

Fairly large relative variations around the average 
sound speed Cg can be achieved experimentally, for in- 
stance by means of a Feshbach resonance that modifies 
the coupling constant g along the flowing BEC. Hence we 
shall consider values of D from 0.1 to 0.7 (higher values 
are not excluded experimentally, but proved difficult to 
reach with our code). The corresponding values of A go 
from 70 to 10, while the order of magnitude of Hawking 
temperature is 



Th ^ 0.1 - 0.8 nK. 



(75) 



^ We are grateful to Eric Cornell for providing us with realistic 
experimental values for all parameters and to Jeff Steinhauer for 
further comments. 
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In general there will be a variation of both the sound 
speed and the flow velocity. Which one varies most de- 
pends on the experimental conditions. In [14], q was 
taken to be zero. In the following, the typical value of q 
will be 0.3. In Sec. IV C 31 we shall nevertheless explore 
the whole range < q < 1. 

In realistic conditions a condensate has an effective 
temperature approximately fixed by the chemical poten- 
tial fjL ~ mcQ. To be able to quantitatively compare 
Hawking temperature with this effect, we define 



UICq 



30 nK, 



(76) 



and in the following we consider initial temperatures (en- 
tering Eq. ([is]) ) 



?0J 



(77) 



with T ranging from 1/3 to 1, so that Tin goes from 10 nK 
to 30 nK. Equation ^ together with ^ yields 



(78) 



If A is constrained by Eq. ([75]) . this can be rewritten as 



T 



30 



D' 



(79) 



The fact that Tin is about two orders of magnitude higher 
than the Hawking temperature will surely complicate 
measuring the spectral properties of the Hawking radia- 
tion. 

In the following sections, we start by studying these 
spectral properties assuming the condensate has zero 
temperature. The analysis when there is an initial tem- 
perature is then performed in Sec. IV D] Our main goal is 
to understand how the adimensional parameters D, A, q 
affect the fluxes. To this end it proves very convenient 
to work with the rescaled frequency uj/k and rescaled 
energy flux defined below. 



1. Typical spectra 

In Fig. [7]/^ and T^/Th are represented versus uj/k for 
D — 0.1, 0.4 and 0.7. The corresponding A are fixed by 
Eq. JTSD and are respectively A = 70, 18 and 10. The 
parameter q specifying the relative contribution of c and 
V to the gradient at the horizon, see Eq. ([H]), is fixed to 
0.3. The three values of Wmax/'* are respectively 1.17, 2.2 
and 2.54. The energy flux and T^^ quickly drop to zero 
when approaching Wmax, as expected. Until short before 
T^ is nearly constant. 

Note nevertheless that this constant temperature can 
differ from Th'- for D = 0.1, the asymptotic temperature 
when w ^ is equal to Tq = 1.08 Tj^. For D = 0.4 and 
D = 0.7, we found To = Om&TH and Tq = 1.0004Tff 
respectively, so To differs from Th only by a fraction of a 
percent when D is large enough. Note also that the scale 
separation condition A ^ 1 is not the relevant criterion 
to predict the importance of the deviation with respect to 
the standard spectrum, since A is much larger for T) = 0.1 
than for the other two values. Instead, it is the ratio 
(^max/^ that controls the deviation from the standard 
temperature. This is conflrmed in Sec. IV C 21 

The thermality of the spectra can be characterized 
more precisely. To this end we define 



Tn-T 



LU = kBTo/h 



To 



(82) 



that measures the running of the temperature. For the 
three spectra of Fig. [71 we found |AtoI — 0.1%. This 
establishes that, to a very good approximation, the en- 
ergy spectra are Planckian, but truncated very close to 
(^max- At this point it should be stressed that to reach 
this result, nowhere have we used the gravitational anal- 
ogy. In fact the analogy is validated by our results since 
it would have predicted a thermal spectrum at the stan- 
dard Hawking temperature, missing however the trunca- 
tion near Wmax and the small running IAj-qI). 



C. Spectral properties of HR at zero temperature 

Let us first study the energy flux of positive frequency 
u quanta on the right of the horizon, i.e., what corre- 
sponds to Hawking radiation. We denote by F the total 
energy flux, and we deflne the energy flux density as: 



2tt dF „ |2 
kb-Lh duj K 



(80) 



The factor 2Tr/ ksTH is here for convenience, so that f^, 
be dimensionless and normalized to 1 at cj = when the 
occupation number n^^ — is the standard, Planckian 
one with temperature Th. To characterize the deviations 
from the standard flux, it is also convenient to use the 
effective temperature T^ deflned as 



1 



exp{huj/kBT^) - 1' 



(81) 



2. oJmax/K controls the deviation from the standard 
spectrum 

In this section, in order to show that Wmax/'* is the 
main quantity that controls the modiflcations with re- 
spect to the standard Planckian spectrum with temper- 
ature Th, we relax the constraint (|73p . and allow for 
arbitrary values of A for a given D. We characterize the 
deviation with respect to the standard result by the rel- 
ative difference 



J UJ J to 



fH 

J UJ 



(83) 



as a function of cj,„ax/'*, where denotes the Planck- 
ian energy flux density with temperature Th and where 
LUH = kBTn/h. 
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Figure 7: Energy flux density fui (left plot) and elTective temperature T^i/Th (right plot) versus ui/n, for {D,X) = (0.1,70), 
(0.4, 18) and (0.7, 10). q is flxed to 0.3. 




Figure 8: Ah versus Wmax/ft: for D = 0.1 and D — 0.7. q is 
flxed to 0.3. 



In Fig. [H for q = 0.3, the curves associated with 
D — 0.1 and D = 0.7 have similar shapes, with a maxi- 
mum around lo^i^^/ k = 0.6, with height 0.35 for D = 0.1 
and 0.18 for D = 0.7. The deviation then decreases and 
in both cases becomes less than a percent as soon as 
^max/^ > 2. This confirms what was found (in a differ- 
ent setup) in Rcf. [17]. This is not the end of the story 
however, since as we now show, there can be significant 
deviations from thermality in a BEC, even when Wniax/'* 
is large, depending on the value of the parameter q. 



3. Effect of inhomogeneity: deviations from thermality 

In the preceding sections, we fixed q — 0.3, thus re- 
stricting attention to experimental cases where the varia- 
tion of c+v is mainly due to the sound speed. This would 
be the case for instance if one used a Feshbach resonance 
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Figure 9: Effective temperature as a function of uj/n for D = 
1, q — 0.7 and A = 6. This situation describes approximately 
the experimental realization of Ref. [25l ]. 



to vary the coupling constant g across the BEC. In [25|, 
the authors report to have created a WH/BH horizon pair 
in a BEC using a different technique, where a local in- 
crease in the flow velocity leads to a decrease in the speed 
of sound. Their profile c+v is not symmetric with respect 
to the horizon, contrary to our parameterization, and has 
two horizons. Ignoring the WH horizon we can approxi- 
mately describe their experimental realization within our 
setup. From their experimental values one gets A ~ 6, 
Z? ~ 1, and q ~ 0.7, which yield (j-'max/'* — 3. Given 
the results of the preceding section, one expects a robust 
spectrum of phonon radiation. This is however not the 
case, as shown in Fig. [9] where the effective temperature 
is shown as a function of the frequency. The temper- 
ature becomes equal to Th only for frequencies w ~ k. 
In the low frequency part the power is suppressed. From 
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Figure 10; Effective temperature as a function oiu/ n for D = 
0.4 and (g, A) = (0, 19.54), (0.5, 17.15), (1, 15.47). tJmax/«: is 
the same for all curves, equal to 2.2. 
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Figure 11: \A^\'^ (upper plot) and versus lo/k. The pa- 
rameters are identical to those of Fig. 1101 



this case we learn that the precise mixture of c and v used 
to form the horizon significantly affects the properties of 
the Hawking flux and that the radiation produced in the 
setup of Ref. [1^ should not have a thermal spectrum. 

The effect of q is investigated more systematically in 
Fig.Hni for a lower value of Wmax/^ so as to ensure better 
numerical control. In that figure, T^/Tu is plotted versus 
lo/k, for D = 0.4, and q — (homogeneous BEC where 
only c varies), q — 0.5 and q = 1 (only v varies). For each 
value of g, A is tuned so that we work at fixed LOmax/ 
to facilitate the comparison. For g = and q = 1, the 
effective temperature varies significantly and the radia- 
tion is thus not thermal. For q — 0.5, as in the previous 
plots with q — 0.3, the radiation is (almost) thermal. We 
verified that the same results hold with other values of 
D and A, with the same particular role of q = 0.5. 

The strong effect of q on the Hawking fiux can be un- 
derstood as follows. Contrary to the setup studied in [l3| 
where the u and v sectors were completely decoupled in 
the dispersionless limit, there is no such decoupling in 
the present settings. The creation of left-moving quanta 
and the elastic scattering of u quanta into v quanta will 
thus be large, even when the dispersion plays no role. 
The above results thus suggest that for q close to 0.5, the 
u-v mixing is small and gets larger for extreme values 
of q near and 1. To confirm this, we now turn to the 
properties of the flux of left-moving quanta, fi^ . 

4-. Mixing of u and v quanta 



at W = Wmax- 

Near Wmax, goes to zero in the three cases. Instead, 
the low frequency behavior of changes dramatically 
depending on the value of q. First, the curves correspond- 
ing to g = and q — 1 join at low frequencies, and are 
proportional there to \/lo. This behavior in 1/uj means 
that the energy flux Hlj carried by the v quanta is con- 
stant and nonvanishing at low frequencies. Secondly, for 
q — 0.5, stays everywhere below 10"'^ and is propor- 
tional to Lo at low frequencies. This proportionality with 
uj was also obtained in the different setup of Ref. pTf . 
and seems to indicate that the case q = 0.5 (and not 
q = 1 as one would have thought naively) is effectively 
similar to the setup of that reference. Even though this 
particular point is most probably an artefact of Eq. 
where v and c follow the same function, the qualitative 
conclusion that the u-v mixing is lower when the horizon 
is formed by a variation of both v and c rather than of 
only one of them, should hold experimentally. We em- 
phasize this point because a small u-v mixing guarantees 
that the Hawking B type of correlations is the largest, as 
explained at the end of Sec. IIV F 21 

Figure [T^] shows the influence of D on the u-v mixing 
coefficients and for q = 0.3. Both the creation 

of v quanta and the elastic scattering coefficient A^^ are 
significantly affected by D and increase by more than 
one order of magnitude between D — 0.1 and D — 0.7. 
This will have important consequences when taking into 
account an initial temperature. 



Figure [TT] shows the particle flux of v quanta and 
the elastic scattering coefficient |v4„P as a function of 
uj/k for the same parameters as in Fig. [TOl 

At low frequencies, \A^\'^ is nearly constant. It is two 
orders of magnitude smaller for q = 0.5 than for the other 
two values. It should also be noticed that \A^\'^ does not 
vanish for lu Wmax since both 0" and 4>^ remain well- 
defined above usniax so that \A^\'^ connects smoothly to 



5. Dispersionless elastic scattering 

It is instructive to look at the behavior of lA^^^jp eval- 
uated for hiOH = ksTn as a function of Wmax/'*, and for 
fixed values of D and q. It is shown in Fig. [T31 When 
^max/ K is greater than about 2, becomes nearly 

constant, with a value that depends on D and q. This 
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Figure 12: \A^\'^ (upper plot) and versus uj/k, for various 
values of D. The legend applies to both plots. The values of 
A are identical to those of Fig. [7] and q is fixed to 0.3. 




Figure 13: evaluated for huj = ksTH, versus ujrna,x/n for 

D = 0.1, D = 0.4 and D = 0.7. g is fixed to 0.3 in the upper 
plot and 0.5 in the lower one. 



means that for large Wniax/^, the dispersion no longer 
plays any role. This is to be opposed to what was found 
in [l7|, where jA^^^P scaled approximately as A"''. The 
reason is that in a BEC, the wave equation Eq. (|16p does 
not factorize into a u and a v part in the dispersionless 
limit A —^ oo. Thus, there is always some elastic scatter- 
ing between u and v modes. What Fig. [01 shows is that 
the dispersionless value of l^cuff P is quickly reached, or 
equivalently that the A-dependent contributions vanish 
rapidly for increasing A. This A-independent elastic scat- 
tering exists also for q = 0.5, whereas we have seen that 
in this cas e n^, has a behavior similar to the correspond- 
ing one in [17|]. This is somewhat surprising. 

Note finally that, since \A^\'^ is nearly constant in cu, 
see Fig. [TT] and Fig. [T^l the value for some {D,q) 
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Figure 14: Effective temperature divided by Hawking tem- 
perature as a function of lj/k for q = 0.3 and {D,X) = 
(0.2, 3.5) (solid line) and (0.4, 1.8) (dashed line) when the vari- 
ation oi c + V occurs on one healing length. 



in the regime where A plays no role actually gives l^^jp 
for all oj. 



6. Strongly dispersive regimes 

In Fig. in we saw that, starting in the robust regime 
(^max/K > 2 and reducing Wmax/'*, first increases, 
which indicates also an increasing asymptotic tempera- 
ture To J E^nd then monotonically decreases for uj^i^^/ k < 
0.6. We refer to the latter regime as the strongly dis- 
persive regime. It is reached when the variation of c -f w 
occurs on much smaller distances than what is assumed 
in Eq. dZl. 

In Fig. [lH the variation is assumed to occur on one 
healing length. This amounts to change the numeri- 
cal factor in Eq. dTS]) so that now A/k = 0.7/D. The 
effective temperature Ti^ is shown for q — 0.3 and 
{D,X) = (0.2,3.5) and (0.4,1.8). The corresponding val- 
ues of Wmax/«^ are 0.16 and 0.22, that is, much smaller 
than in Fig. [T] The shape of the two curves is remarkably 
similar to what was obtained in the right plot of Fig. [71 as 
soon as the frequency is less than about 0.2 x o^max, is 
nearly constant (the running A^^ is equal to 7% and 5%) 
and then drops quickly to zero when approaching Wmax- 
However, Tq, the asymptotic value of T^, significantly 
differs from the Hawking temperature: it is respectively 
equal to 0.6 Th and 0.75 Th- From this we conclude that 
when dispersion is strong, the low-frequency spectrum 
remains in 1/w as for a Planck spectrum. However since 
the cutoff frequency Wmax is only ftwmax/l^BTo) = 1.7 
and 1.9 respectively, the spectrum is no longer Planck- 
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Figure 15: Effective temperature T^^ divided by the healing 
temperature of Eq. ((76}, as a function of oi/wmax for q = 0.7 
and D = 1 in extremely dispersive cases, when k ^ cq/^o- 



ian. ^ We finally notice that T^/Th becomes smaller and 
smaller as ujraaxi decreases. 

It is thus interesting to consider extremely dispersive 
cases where CqK"^ is much smaller than ^q- In Fig.fTSlwe 
represented T^^/T^^ as a function of w/wmax, with D = 1, 
q = 0.3 and small values of A from 0.1 to 10"'^, corre- 
sponding to w,nax/K ranging from 4 x 10^^ to 4 x 10^*. 

still tends to a constant value Tq at low frequencies 
and the low- frequency part of the spectrum still behaves 
as uj~^. The ratio fiwmax/ {ksTo) is nearly constant, equal 
to 1.6 for all three curves. The running is large, of the or- 
der of 12%. An important result is that the ratio Tq/T^^ 
saturates at a constant value as Wmax decreases. ^ This 
demonstrates that k becomes irrelevant and that Tq is 
fixed only by ^o- We have verified that this result also 
holds for other values of D and q, but with different 
asymptotic values for Tq / T^^ . 



D. Spectral properties with an initial temperature 



Energy spectrum 



As pointed out in Sees. IIVBI and IV Bl the residual 
temperature of a condensate is expected to be about 
two orders of magnitude higher than Th- The spectra 



^ We thank T. Jacobson for a remark about this point. 

^ This result is corroborated by [48] where a step-like variation (in 
a;) of the sound speed is considered in a homogeneous condensate 
flowing with a constant v so as to get a transition from a sub- 
to a supersonic flow. This case can be handled analytically and, 
in our language, it corresponds to the limit k — » oo, with D, 
flxed and q = 0. 



of the previous sections are thus unlikely to be observed 
as such and it is necessary to include the effects of an ini- 
tial temperature. This is easily done using Eq. (j47| and 
the numerical values of the coefficients of the Bogoliubov 
transformation. Assuming that the three initial occupa- 
tion numbers are characterized by a common comoving 
temperature, they are given by Eq. (|45p . Their calcula- 
tion reduces to the computation of the functions ri*"(a;) 
for the three types of modes, which is easily done by solv- 
ing Eq. Since SI™'" are three nontrivial functions 
of Lu, the initial distributions are not Planckian in uj. 

The energy spectrum emitted to the right of the hori- 
zon, with a nonzero initial temperature, is defined as 



fin 



2n—n 

K 



fin 



(84) 



It is represented in Fig. \W\ for the same set of parame- 
ters as in Fig. [7l Two values of the initial temperature 
are considered: a conservative one. Tin = 30 nK, and an 
optimistic one. Tin = 10 nK. With the notations and 
the constraints of Sec. IV B[ they correspond to r = 1 
and T = 0.3 respectively. For D = 0.1, D = 0.4 and 
D = 0.7 the ratio T.n/Tn is respectively 300, 75 and 43 
when r = 1, and 90, 23 and 13, when r = 0.3. 

The spectra differ greatly from those at zero tempera- 
ture, and have a nontrivial behavior. Without surprise, 
they no longer vanish when approaching Wmax because 
neither ri*" nor \Ai^\'^ do. The most interesting point 
is that for D = 0.4 and 0.7 and the lower value of r, 
the spectra follow relatively closely the zero-temperature 
ones until w ~ k. Thus, for frequencies below k, large 
values of D, and low initial temperatures, the measure 
of the phonon energy spectrum to the right gives a good 
estimate of the zero-temperature flux (i.e.. Hawking ra- 
diation) . 

To better understand these spectra, we show in Fig.fTTl 
the contributions to the energy flux of each of the three 
terms in Eq. (|T7l) . along with the full spectrum, in the 
two extreme cases D — 0.1 with t = 1, and D = 0.7 with 
T = 0.3. 

In both cases, at high frequencies, the initial distribu- 
tion largely dominates the spectrum, and all we see 
is just the flux of the initial quanta. At low frequencies 
instead, the main contribution comes from the sponta- 
neous plus stimulated emission term |/3(^p(l-|-n™-|-n™^), 
cafled "creation term" in Fig. [iTl For D = 0.1, the 
u — V scattering term is small and never dominates, but 
it is never completely negligible. The contribution from 
l^wP^i™'" quickly becomes comparable to the exponen- 
tially decreasing stimulated emission. 

For D = 0.7, this contribution is not much smaller 
than the creation term at low frequencies, as could be 
expected from the results of Fig. [T^l On the other hand, 
the contribution from n™ long remains very small com- 
pared to the other two, until about uo = k. Thus the 
creation-|-stimulated emission term is actually dominated 
by the Hawking pair creation effect, since fi™^ is close to 
n™. This explains the relative similarity between the 
spectra with r = 0.3 and r = 1 in Fig. [111 when D = 0.7. 
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Figure 16: Energy flux emitted to the right of the horizon, q is flxed to 0.3. The dotted spectra are those of Fig. [T] The values 
T = 1 (solid lines) and r = 0.3 (dashed lines) correspond respectively to an initial temperature Tin ~ 30 nK and Tin ~ 10 nK, 
see Eq. (177)1 . 
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Figure 17: Different contributions to the energy flux to the 
right of the horizon. Left plot: D = 0.1, r = 1. Right plot: 
D = 0.7, r = 0.3. 



In brief, the main lesson from these plots is that, in 
general, there is no clear hierarchy between the various 
contributions in n^"^, which makes the interpretation of a 
measurement highly nontrivial. Nevertheless, let us try 
to identify what the optimal conditions are for such a 
detection. 



2. Optimal experimental conditions 

The comoving frequency Q of the right-moving quanta 
is strongly redshifted in the black hole geometry, which 
implies that the initial f7"'™((jj) is much larger than uj. 
Since the integrated rcdshift increases with A, one could 
ask whether lowering the Hawking temperature with re- 
spect to the initial temperature Ti„, with a fixed value 
of D, could not improve the results above. Indeed, if the 
increase of the redshift were such that W^"'/Tin would 
grow with A, the initial occupation numbers n™ and 
ri™^ would decrease exponentially, and one could hope 
to have the stimulated plus creation term dominate over 
the other two on a larger interval of frequency uj. This is 
not the case however, as we now explain. 

The initial proper frequency n™{uj)/K for the three 
types of modes is shown in Fig. [THl for D = 0.4, q = 0.3, 
and different values of A. As expected, at a given uj, 
is an increasing function of A. However it scales 
only as A, and since the ratio Tm/Tn also is proportional 
to A, n^'^''^/Tin hardly changes. On the other hand, the 
left movers suffer almost no redshift, and O™'" has al- 
most no dependence on A. Thus, D,™''" /Tin decreases as 
A^^ and the initial occupation number of the v modes 
increases as A for a given u. 

These remarks are summarized in Fig. I19[ where the 
occupation numbers n™ of Eq. (|45p are represented for 
the same parameters as in Fig. [THl and fi™^ are al- 
most equal and, more importantly, do not change with 
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Figure 18: Initial proper frequency fl^" /k as a function of uj/k for the u, u> (left plot), u, —uj (middle plot) and v, u> modes (right 
plot), for D = 0.4, q = 0.3 and t = 0.3, and different values of A. The legend applies to all plots. 
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Figure 19: Initial occupation numbers n^" (left plot), n'^^ (middle plot) and n^"'" (right plot) versus oj/k, for the same values 
of the parameters as in Fig. [18] 



A at low frequencies, n™'" on the other hand explodes 
when A increases. Since \Auj\'^ does not depend on A 
in first approximation, the consequence is that the scat- 
tering contribution in n£'" dominates at low frequencies 
when A increases, since the other contributions remain 
unchanged. In conclusion, there is no gain in lowering 
Th- 



In fact, as can be seen from the figures of the previous 
section, the optimal conditions to be able to measure 
directly the zero temperature flux fuj are reached when 
there is a large relative variation of c + w on the smallest 
possible distance, that is the largest possible D with the 
smallest possible A. Remember also that c and v should 
contribute nearly equally to the variation for to be 

as small as possible. 



3. Density fluctuations and atom flux 

As explained in Sec. IIV C[ the flux /^^ studied in the 
previous sections is observable only if one can distinguish 
experimentally between left- and right-moving phonons. 
It is thus important to address the question whether the 
easily accessible density fluctuations also contain some 
signature of the particle creation process. The function 
PqvG^^ of Eq. ((52| . evaluated in the coincidence point 
limit in the right asymptotic region, is represented in the 
upper plot of Fig. [101 for D = 0.7, q = 0.3 and r = 0.3. 
This quantity is adimensional, does not depend explicitly 
on pq and cq and is proportional to the power spectrum of 
the density fluctuations. In the lower plot, .F™™, deflned 
in Eq. (|56p and adimensionalized by multiplication with 
Povcq/k, is represented, for the same set of parameters. 
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Figure 20: Upper plot: adimensionalized (by multiplication 
by pov) power spectrum of density fluctuations as a function 
of ui/n. Lower plot: adimensionalized comoving atom flux as 
a function of u/k. In both plots, D = 0.7, q = 0.3 and r — 
0.3, and GJ^ and are evaluated in the right asymptotic 

region, in the coincidence point limit. "HR related" refers to 
the contribution weighted by n£'" and "initial excitations" to 
the one weighted by nZ'^" ■ 



In both cases, the contributions governed respectively by 
n^" (referred to as "HR related" in the figure) and by 
n™'" have been plotted separately. They are of the same 
order of magnitude but for cu <ti k, that is, in the re- 
gion where n^" is mainly due to Hawking radiation, see 
Fig. [ini the contribution from the initial excitations is 
greater by a factor of about 2 in the density fluctuations, 
and of up to 6 in the atom flux. This confirms the remark 
made in the theoretical analysis that no clear signature of 
Hawking radiation can be observed in these local observ- 
ables. However, as pointed out in Sec. IIV E[ one could 
hope to combine both observables to extract n^"-. 



VI. CONCLUSIONS 

In this work, we presented a complete description of 
the scattering of the phonon modes propagating in sta- 
tionary, one dimensional condensates that possess a sonic 
horizon. Our description is based on the BdG equa- 
tion, and the scattering in the horizon region is expressed 
in terms of 3 x 3 Bogoliubov transformations relating 
asymptotic one phonon modes. 

We have shown that this scattering affects two types of 
observables, local ones such as the density fluctuations, 
and the long distance correlations. The former are gov- 
erned by expectation values which are diagonal in the 
occupation number, see Eq. (|52|) . whereas the latter are 
determined by interference terms, Eq. (j59p . which reveal 
the entangled nature of the final state. 

When taking into account the condensate temperature, 



the local observables are in general dominated by the ini- 
tial distributions, whereas nonlocal correlations are am- 
plified by these initial distributions without having their 
pattern modified (when the initial state is uncorrelated) . 
Therefore the latter probably constitute the clearest in- 
dication that the Hawking effect is taking place. How- 
ever, it is worth noting that the spatial structure of the 
pattern is not specific to HR strictly speaking, i.e. to a 
nearly Planckian spectrum, since it is determined before 
all by the structure of the 3x3 transformation of Eq. ([1T|) 
rather than by the value of its coefficients. In fact, simi- 
lar patterns are found when considering the scattering of 
classical waves, as explained in App. [Cl and in the lim- 
iting cases At — > oo and k = (with nonvanishing higher 
derivatives of c -|- at a; = for the latter) where the 
spectrum is no longer approximately Planckian (as we 
verified) . 

In the last part, we numerically integrated the BdG 
equation, and obtained the occupation numbers of the 
three kinds of phonons, both with and without an initial 
temperature. The main results are the following. 

Firstly, in the initial vacuum and when the u-v mix- 
ing is small {q ~ 0.5), the spectrum of the outgoing 
(right-moving) phonons is Planckian, with a tempera- 
ture determined by the gradient k of Eq. ((20)) . as soon as 
i^max/K ^ 2 and for frequencies u < o^max, see Fig. [T] 
This result, derived directly from the BdG equation 
(without further approximation), makes precise the do- 
main of validity of the analogy between relativistic fields 
in black hole metrics and a phonon field in the corre- 
sponding "dumb hole". When a^max/^ < 2, the analogy 
fails even for the lowest frequencies, see Fig. [141 which 
invalidates the condition usually found in the literature 
that the analogy is applicable for wavelengths greater 
than the healing length. When the u-v mixing is not neg- 
ligible, the spectrum deviates from the Planck spectrum, 
see Figs.fTUlandfm This is due to grey body factors (that 
can also be computed to a good approximation using the 
gravitational analogy [1^). 

Secondly, as for the detection of the analog HR in 
a BEG, if one can distinguish left- from right-moving 
phonons, the Hawking flux could be observable when 
the relative variation of c + u is large and happens on a 
small number of healing lengths. Instead, if they cannot 
be distinguished, local observables are completely domi- 
nated by the initial temperature of the condensate. One 
must then resort to nonlocal observables such as density- 
density correlations [H, [13] . 

Thirdly, in App. |D1 we established the relationship be- 
tween the fluxes emitted by white and black hole flows. 
Because of the blueshift effect in WH flows, the correla- 
tion pattern is more amplified by the presence of initial 
quanta than in BH flows and should therefore be easier to 
detect. This could well be the optimal case for observing 
HR in the laboratory through the correlation pattern. 
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Appendix A: NONSTATIONARY CASE 

In this appendix we consider nonstationary conden- 
sates. As we shall see, little is modified when compared 
with the stationary case studied in the body of the pa- 
per. To ease the comparison we still work with the one 
dimensional case. The extention to the three dimensional 
case is trivial. 

In the mean field approximation, the condensate is now 
described by 



(Al) 



This wave satifies Eq. (O where V, and g depend both 
on X and t. (As pointed out in [lj|, some tuning of V 
and g might be necessary not to produce time-dependent 
effects that might hide the Hawking radiation.) The con- 
servation of the number of atoms follows from Eq. ([5]), 
and gives the continuity equation 



dtPa + dxifov) = 0, 



(A2) 



where v(t, x) = hko{t, x)/m is the velocity of the conden- 
sate, and ko{t,x) — dxWo{t,x) its wave vector. Plug- 
ging Eq. ((XT|) into Eq. © and using Eq. ((Ml) still 
gives Eq. ([5]) where the chemical potential is replaced 
by hu!o{t,x) — —hdtWo. Then, because of Eq. (|A2|) . the 
condensate is still characterized by v{t, x) and the speed 
of sound 



c^{t,x) = g{t,x) 



Po{t,x) 



(A3) 



To describe the phonon modes, as in Sec. |TT1 it is con- 
venient to work with the relative fluctuation (f> defined 
in Eq. (HT))- Then, using Eqs. jJTl ESI, one still gets 
Eq. p2| . To get the c- number modes which shall be 
used to proceed to the second quantization, we look for a 
complete orthonormal family of doublets Wj — {(j)j,ipj), 
where the index j now does not correspond to a con- 
served quantum number. The scalar product with re- 
spect to which orthonormality is defined is still given by 
Eq. (|77)) . with po{x) replaced by po{x,to). (The choice 
of the time to when the product is evaluated does not 



change its value since it is conserved.) This family enters 
(j) as in Eq. (flT 



(f>it,x) 



E 



4'jit,x) + a] {ipjit,x)) 



(A4) 



where aj and a| are annihilation and destruction oper- 
ators, with which one can construct the phonon Fock 
space. However, contrary to the case studied in the text, 
there is in general no clear interpretation of the vectors 
in this Fock space in terms of particle content, and the 
notion of phonons is inherently ambiguous. The notion 
of phonon can be recovered if the condensate is station- 
ary in the asymptotic past and the asymptotic future. In 
that case, as in the text, one can define two Fock spaces, 
the in one and the out one. 

Inserting Eq. (jA4p in Eq. p2)) and taking the commu- 
tator with dj and a| yields: 

[ih{dt + vdx) - Tp - mc^] (jjj ^ mc^ ipj, 
\-~ih (dt + ivdx) — Tp — mc I (pj = mc (f)j. 

As in the body of the paper, one can eliminate (pj and 
obtain 



[n [dt + vdx) - iTp] 1 [h {dt + vdx) + iTp] 



'2mT, 



(A6) 



= 0. 



We notice that all kinetic terms are of the form of 
Eq. p^ . This is related to the fact that Tp is self-adjoint 
when using the scalar product of Eq. (P7)) . 



Appendix B: ADDITIONAL REMARKS 
CONCERNING EQ. ([16)) 

1. Eikonal approximation and importance of 
ordering 

It is worth exploring the properties of Eq. (fTB)) . the 
stationary version of Eq. (|A6p . We first notice that 



1 



[h {lj + iv dx) + Tp] — = [h {uj + iv dx) + Tp 



{ihv dx +Tp) ,— 



(Bl) 



This makes explicit that the spatial gradient of affects 
Eq. (fTO)) only through a commutator. We also notice that 
•p^, the other mode of Eq. (fT4)) obeys 



I [-h (uj + iv dx) + Tp] -^[h{uj + iv dx) + 

- h^vdx-dx\pu> = 
V J 



To 



0. 



(B2) 



The only differences between the equations for (jji^ and 
for come from the sign of the commutator between 
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V dx and c~^. When this commutator is neglected, as it 
is in the WKB approximation, (/)^ and (p^ thus obey the 
same equation. 

When working to leading order in a WKB (or eikonal) 
approximation, i.e., inserting ~ ^ijdxk^ ([TB|) 
(or inserting (p^^ ~ ^ijdxk^ j^-^ (|B2[) ) gives the disper- 
sion relation in a moving fluid of velocity w. 



2„2 



(cj - kvY = VL^ = k^c 



Am? 



= c^e 1 1 



(B3) 

where ^ = h/^/2mc is the healing length. We recover 
the quartic dispersion relation between the frequency in 
the comoving frame = (w — fcu), and the wave vector 
k. It should be noticed that in a stationary nonhomo- 
geneous flow, 17 and k depend on x through v{x) and 
c(x), whereas the frequency w is a globally defined con- 
stant. Remember also that oj is not necessarily the lab 
frequency, because it is defined in the frame in which the 
condensate quantities only depend on x. 

Returning to Eq. (|16p. it is clear that its role is to fix 
the exact properties of the ODE obeyed by (j)^;. These 
properties could not have been inferred starting from 
Eq. (|B3|) . and applying the substitution k —id^, be- 
cause this naive rule could neither predict the ordering 
of r, v{x) and c?{x) found in Eq. ^6]) . nor the different 
one found in Eq. (|B2p . 



2. Hydrodynamical limit and Euler mode equation 

In the hydrodynamical (dispersionless) limit, for long 
wavelengths with respect to the healing length ^, one can 
drop the two operators Tp in Eq. Then Eq. ITBl) 

reduces to the Eulerian mode equation which describes 
sound waves in a moving fluid, and this with the correct 
nontrivial ordering of T, v{x) and c?[x). Let us verify 
this. 

The Eulerian action is usually written in terms of the 
velocity potential ?/;, related to the velocity fluctuation 
by 5v = dxip, see e.g. p3|. To make contact between ^ 
and the (j) held, one should compare the fluctuations of 

described as in Eq. (|lip to those written as 



«..(../v»>*.(i^ + 

Using dv 



i6W 



(B4) 



^dxSW = dxip, one obtains 
6p 



2i 



Po 



= = SW 

n 



(B5) 



Using the phase fluctuation 6 = rnip/h rather than -(/; 
itself, the Euler action is 

Se^^ I dtdxpo m + vdx) ef ~ {dxBf 
Irti J yc^ 

'(B6) 



Then, using the continuity equation, the mode equation 
reads 



1 



1 



[dt + vdx) {dt + vdx) dx Po d.. 



= . (B7) 



When working at fixed w = idt, and using vpo = est., 
as announced, one recovers the dispersionless limit of 
Eq. (fTB)) obtained by taking the limit Tp ^ 0. 

Notice that Eq. (|B7p is a generalization of the (disper- 
sionless limit of the) mode equation which has been gen- 
erally studied in the literature, see Refs. [sl. H. [sl. [itI. . 
In those references, the equation also followed from 
Eq. (|B6p . but with the extra hypothesis that both c and 
Pa can be approximated by constants, in which case one 
recovers the massless relativistic 2D mode equation. 



3. Link with Unruh's dispersive models 

It is worth noticing that the dispersive properties of the 
phonons come through the two operators Tp in Eq. (fTB|) . 
This is not what we would have obtained had we used 
the rules of [3| with a quartic superluminal dispersion. 
Indeed, the Eulerian action supplemented by a quartic 
term {k^ / K^) is 

5a = — / dtdxpoi^ [{dt + vdx) ef 



{dxef + Y^{dlof}^ 



(B8) 



and the corresponding mode equation reads, in a station- 
ary condensate and at fixed w. 



{uj -f ivdx) -^{uj + ivdx) - V dx- dx 



-d'-d' 

A 9 ^X ^X 



V 

^ 0. 



(B9) 



In nonhomogeneous situations, the quartic term encod- 
ing the dispersion differs from that of Eq. (|16p . which 
means that 0,^ [x) will differ from (j),^ [x) in any nontrivial 
background. 

Finally, it is also worth noticing that Eq. can be 
obtained from an action for a single field with a quadratic 
kinetic term, and which generalizes the Euler action: 



Sd, — 



— / dtdxpa{ i 



T 

vdx + 

n 



[dx^f 



(BIO) 



This action can be obtained from S = 
J dtd^x(ifi^*dt'^ - H), where H is given in Eq. (O, 
using Eq. (fTT|) , keeping all quadratic terms in cj), </>* , and 
using Eq. (fT^ to eliminate (j>* in favor of (j>. The conju- 
gate momentum that enters the equal time commutator 
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(x),7r(y)] = ihd{x — y) is rather unusual 



ihpo [(j) + 0^) = ihSp, 



where we used Eq. (jB5|) . It thus obeys vr 
Eq. (|A4p . TT decomposes as 



(Bll) 
-tt'''. Using 
■ (B12) 



where tTj = ihpo{(j)j + ^j) — ^T^j- The conserved scalar 
product (which generalizes the standard Klein-Gordon 
one [1^ and which agrees with Eq. ((27| ) is 



(B13) 



It is not clear whether this alternative way of describing 
the phonon field presents any advantage over the origi- 
nal version of Eq. (fT5|) . It could nevertheless be useful 
to study how dispersion affects the analogy with gravita- 
tional systems. [H, H^, 



Appendix C: CORRELATION PATTERNS IN THE 
CLASSICAL LIMIT 



Rather than using the vacuum or a thermal state as 
in the text, we assume that the initial state also con- 
tains a highly excited coherent state.'' This state can be 
obtained by making use of the displacement operator 



exp [w 



(CI) 



Any of the three initial operators a"'^ appearing in 
Eq. (j5ip can be used to construct the corresponding ini- 
tial wave. One can also consider wave packets engendered 
by 



exp 



dhj {w al - w* a^) 



(C2) 



When imposing the normalization J duj\d^\'^ = 1, the 
mean occupation number of initial quanta added by Dd^ w 
is jwp, as it is for in Eq. (jCip . 

When p™, the initial state without the coherent state, 
is such that Tr[/3"^ x(^i 2:)] = 0, the anticommutator in 
the presence of the coherent state separates as 

1 



Git,x;t',x') = -Tr 



Dd,^urDl^') {X{t,x),x{t',x')} 
= Xdit,x)xd{t,x)+G''^{t,x;t',x'), (C3) 



We are grateful to W. Unruh for bringing our attention to co- 
herent states while one of us, R.P., was presenting this work at 
the Peyresq conference of Cosmology in June 2009. 



where G™ is given in Eq. (|49|) , and where the mean value 
of the field operator is 

Xdit,^) = ^^[{Dd,n.rDl^)x{t,x) 

= {0^n\Dl^Xit,x)Da,^\0^n), (C4) 

as if the initial state was the in vacuum. (To get these 
equations, we have used Tr[p™ x(^j 2;)] = and the re- 
lation b\ ,^auibd.w = a^j +wd^, see e.g. [HI.) From 
the decomposition Eq. (jC3[) . we see that the correlation 
pattern is the sum of the pattern encoded in the state p™ 
we studied in the former sections, plus the new pattern 
encoded in the real wave packet Xd{t, x) associated with 
the coherent state. 

Assuming that this wave packet is initially made only 
with a™'", it is given by 



Xd{t,x) 



/>OC 

/ du! [wdi^ e~ 
Jo 



*xL"'"(:c)+c.c." 



(C5) 



At early times, before it enters the near horizon region, 
it describes a single wave packet propagating against the 
flow in the region where the flow is supersonic, see Fig. [6] 
As shown in that figure, at late time, it splits into three 
wave-packets. 



Xd{t,x) 



duj \wdi^ 



—iwt out, u / \ 

e a^jXi^ ' (x) 



duj 



wd^ 



■A^x"J''^^x)+c.c. (C6) 



+ ^""'^ dcj {wd^e-^-' [xl^^'^^ix]] * + c.c.} . 

The first one describes the transmitted wave. It is am- 
plified with respect to the initial wave if J duj\ai^dui\'^ > 1 
which needs not be always the case, as can be seen from 
Eq. ([55]) . The second wave describes the left moving 
packet obtained by elastic scattering. The third wave is 
due to stimulated pair creation process, see the minus 
sign in Eq. p9p . It describes the partner's wave, and 
is present only if d^ has non-vanishing components for 
u! < Wmax sincc P-u, = for w > Wmax- Noticc that 
the Fourier component multiplying the mode e*"*x°'l5'" 
is the complex conjugate of (3-^wd^. The fact that d* 
appears guarantees that when replacing d^ by d^je*'^*", 
the partner wave with negative uj stays synchronized, i.e. 
it is shifted by the same lapse of time as the two other 
wave packets with positive w. Notice also that the pres- 
ence of the coefficient implies that the partner wave 
function cannot be localized in a region smaller than c/k 
because (3^^) er'^^l'^ for w ^ k. In this we recover the 
width of the correlations to the partner which is found 
when using the correlation function 13, 35]. 

We can now relate the space-time pattern encoded in 
G"", the second term in Eq. (|C3|) . to that encoded in Xd- 
Since these two terms have a completely different origin, 
one might a 'priori think that the two patterns will be 
radically different. However, this is not the case. 
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Two conditions must be met for the patterns to be ap- 
proximately the same. The first concerns the out modes, 
whereas the second concerns the Bogohubov coefficients 
appearing in Eq. (jC6|) . At a given late time t, i.e. af- 
ter a lapse of time At, the chosen wave-packet Xd has 
been scattered, the pattern encoded in Xd is given by the 
three values of x where the waves constructively inter- 
fere. Working far from the horizon, i.e. nAt ^ 1 (see 
also the discussion in the paragraph after Eq. ([621) 1. and 
using Eq. ([M)) . these are the solutions of 



At 



HJ, a 



(x) = At, 



(C7) 



where ujd = J dujuj\di^\'^ is the mean frequency of the wave 
packet. When the dispersive effects are weak at late time, 
as it is the case for uJd ^ c/^o, see the discussion after 
Eq. ([M)) . the solutions of Eq. (jC7[) are essentially inde- 
pendent of LUd- Therefore, in the large x limit, one finds 
the same pattern of correlations whether one considers 
the two-point function at some given time t as a function 
of x' given x, see Eq. (|63)) . or whether one looks for the 
two partner waves when using a packet which contains a 
branch that arrives at x at time t. 

In the above reasoning, we have neglected the uj de- 
pendence of the Bogoliubov coefficients. If their rel- 
ative phase would vary rapidly, and therefore that of 
— P^i^/a^ as well, Eq. (jC7P could receive significant 
corrections that might depend on the value of w^, thereby 
giving different patterns for different wave packets. This 
is not the case. When using a massless two dimensional 
relativistic field, one finds that the pattern is universal 
as argZcj is truly independent of see e.g. Eq. (3.49) 
in [26|. Upon considering dispersive fields, when both 
^d c/Co and k <C c/^q are satisfied, one finds [3] that 
the late time pattern is unmodified even though the early 
pattern is completely different and depends on both the 
value of A = c/^o and the sub- or superluminal charac- 
ter of dispersion [i^, [43|. These results were obtained 
using saddle point and WKB approximations which are 
both reliable in the present regime. The analysis was 
performed with a scalar dispersive field, but it also ap- 
plies to a phonon field in a BEC because the differences 
between the various dispersive models, see Appendix iBl 
are not relevant in a WKB regime. 

In brief, we saw that the late time pattern obtained 
from the scattering of classical waves is (to leading or- 
der in w/A) the same as that obtained from amplifying 
vacuum (or thermal) fluctuations. This correspondence 
follows from the fact that both patterns are obtained 
from the (common) decomposition of in modes into out 
modes. Given this we can address a question which has 
been raised by several of our colleagues: to what extent 
observing the scattering of classical waves could be con- 
sidered as a signature of the Hawking effect? 

A convincing signature 21| would consist in measuring 
the contribution of the negative frequency modes which 
anomalously contribute to Eq. ([39]) . once sent a wave- 
packet composed only of positive w, as it is the case in 
Eq. (|C6p . The experimental difficulty is different for BH 



and WH flows. In BH flows, one has to be sure that the 
initial packet contains no negative modes otherwise their 
final contribution will be dominated by the initial one, 
thereby preventing a neat measure of (3^^- The values 
of Lo with P-uj ^ are such that \uj\ < ujma.x- Given 
Eq. ([23]) . they correspond to an interval of initial values 
for fc(cj): 



= A;(-t^max) < k{uj) < /Cmax = A:(tJmax), (C8) 



where 



max 



Am? I 



(C9) 



The restriction to modes with positive frequency imposes 



k>ko = fc(0) = fi( 



2 )l/2 



By a numerical analy- 
sis, we found that the ratio (fcmax — ko)/ko is of the order 
of 20%, independently of m (or ^), and weakly varying 
with D and q. Thus in principle, one could selectively ex- 
cite (e.g. by Bragg spectroscopy techniques [I^llll) cer- 
tain values of k in that interval. However in experiments, 
there will be a tension between the limited lifetime of the 
condensate and the narrowness of the packets in k space 
which gives broad packets in real space with long travel- 
ing times. Orders of magnitude taken from [s^l and [2^ 
indicate that there could be one wave packet satisfying 
all desiderata. 

For WH flows, because of the time reversal symmetry, 
the situation is the opposite in that there is no difficulty 
to prepare an initial wave-packet containing only positive 
frequencies, whereas the final positive and negative fre- 
quency packets might be difficult to distinguish as they 
will follow similar trajectories [2]| . 



Appendix D: WHITE HOLES 

So far we only considered condensate flows that possess 
a sonic horizon which is analogous to that of a black 
hole. However sonic horizons which act as a white hole 
are closely related to the black hole ones. Indeed, the 
transformation 



v{x) 
c{x) 



-v{x), 
c{x), 



(Dl) 



transforms the profile studied in the text, Eq. (|2ip. into 
a white hole profile, in which the characteristics focus 
forward in time, see Eq. with k < 0. (This equation 
now applies to the left-moving modes, so that the quanta 
mainly produced are left-moving ones.) 

The equivalent of Eq. (|A6[) in this WH profile reads: 



(D2) 



{ [h {dt - vd,) - iTp] 1 [h {dt - vd.,) + zTp] 



■2mT„U, = 0. 
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q=0.3 




Figure 21: Occupation number of the v quanta for a BH (n^, 
solid line) and of the u quanta for the corresponding WH 



under the transformation (ID1|) {r 



dashed line), as func- 



tions of cj/k for q — 0.3 and D = 0.4 (upper plot) and D = 0.1 
(lower plot). 



WH flow, when working at fixed frequency w, the modes: 



it,x) 



,.oui,BH 



{x) 



c)]% (D4) 



(and similarly for (^^™'WH-j form a complete in basis. 
Notice that the in/out character and the u/v character 
of the modes is interchanged between the BH and the 
WH. Indeed, because of the time-reversal symmetry, the 
phase and group velocities of the modes change sign. The 
Bogoliubov transformation relating the in to the out WH 
modes is thus the inverse of (|4H) . This implies that when 
dealing with a WH flow in the initial vacuum, instead 
of Eq. (UH), the mean occupation numbers of emitted 
quanta are 



nT = 1/3- 



n 



WH 



n. 



WH 



(D5) 



q = 0.3 




^ 0.01 1 1 — I 




CO/ K 

Figure 22: Relative difference j j/?^ ^ - | /l/?^ |^ as a 

function of uo. Same parameters as in Fig. 1211 



Under a time-reversal t —t, Eq. (jD2|) becomes the 
complex conjugate of Eq. (jA6|) . This proves that there is 
a one-to-one mapping given by: 

r : </.fH(t,x) ^ <j>Y''{t,x) = [0f (-i,x)]* . (D3) 

The same relation holds between ip^^ et ■ ^^'^ dou- 
blets Wj^^ and W]^^ related by Eq. jDSl thus have the 
same norm for the scalar product (l?7)) . 

With these remarks, it is manifest that in a stationary 



The third line shows that the total number of produced 
pairs is equal to that for the corresponding BH. However 
since ^ |/3_t^p, the repartition of the quanta with 

positive uj into left and right movers differs. 

Before considering the implication of these relations, 
it is worth noting that only one quantity governs all dif- 
ferences between the occupation numbers and the elastic 
scattering in the WH and BH cases. Indeed, Eq. (jPSP 
(or the third equation in (HH)) shows that 

Kf + \B^f^\P-J' + \B^\', (D6) 
while the conservation of the norm of 0«^out,BH ^j? 

0«an,BH gj^gg. 

|«.,P + |^P-|/3-J' = \a.\' + \A^\'-mH= !)• (D7) 
Thus, 

- = \B^\^ ~ \B^\^ = \A^\' - (D8) 

The relative difference between the occupation numbers 
in the BH and WH cases is small whenever the u-v mixing 
coefficients are small, see for instance Sec. V C 3 in 
In the case of a BEC, however, the u-v mixing coefficients 
are not necessarily small and can even grow as 1 /lu at low 
frequencies, see Sec. IV 

To further analyze the difference between WH and BH 
fluxes, in Fig. [211 we have represented IB^]"^ and \Buj\^ for 
q = 0.3 and two values of {D, A). Contrary to \B^\'^ that 
grows as l-B^P is nearly constant at low frequencies, 
so that their difference grows like 1/uj. The consequence 
of this is that the relative difference (|/3j^p — |/3_a;p)/|/?ajp 
tends to a constant at low frequencies, as is verified in 
Fig. [221 where this difference is shown for the same pa- 
rameters as the previous figure. This constant value can 
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be significant. For D = 0.4 it is slightly greater than 1%; 
for D — 0.1, it is very small, of the order of 0.1%. For 
both values of D, the relative difference becomes impor- 
tant near Wmax, but both |/3t^P and |/3-cjp vanish when 
reaching this frequency. 

From Eq. (|D5|) several important lessons can be drawn. 
Firstly, we established that, when starting with the vac- 
uum state, the fluxes of phonons emitted by a WH flow 
are directly related to those of the corresponding BH flow. 

Secondly, from this correspondance, the WH flows ap- 
pear to be as stable as the corresponding BH flows. It 
should be stressed however that this conclusion is reached 
when using real frequencies oj, and not taking into ac- 
count the modes that are not asymptotically bounded. 
Whether this implies that WH flows are stable is a moot 
point which deserves further study, see (53 | . 

Thirdly, when expressed in the asymptotic regions in 
terms of the comoving frequency n°'(uj) (or the wave 
number k°'{u!)) the spectral properties of the WH fluxes 
are radically different from those of the corresponding 
BH. This is because, when considering a WH flow, the 
relationships between the comoving fi" and the conserved 
frequency lo are time reversed with respect to those for 
a BH. Hence, unlike what was found for a BH flow, 



Eq. (|46)) now determines the flnal values of the 17° which 
are thus blueshifted. In the ideal case of an initial vac- 
uum, this would be a great advantage since the detec- 
tion of higher frequency phonons is easier. When taking 
into account the nonvanishing temperature of the con- 
densate, the initial occupation numbers are now given 
by n'^''">WH(^) ^ l^gO— ■BH(^)/T„, i^gj^^g they 

are much larger than those in the corresponding BH case 
since fji^out.BH ^ j^. jj^^piigg ^j^^^ ^j^g initial occupa- 
tion numbers in WH flows will always be larger than the 
number of quanta spontaneously created. From this, two 
conclusions can be drawn. First, the detection of the 
analog Hawking radiation, based on the measurement of 
final occupation numbers, will be more difficult than in 
the corresponding BH. Second, on the contrary, the ob- 
servation of the long-distance correlations of Eqs. fST]) 
and (|58p will be greatly facilitated as the larger initial 
occupation numbers increase their amplitude. 

Fourthly, WH might also be more appropriate than 
BH to manipulate coherent states because at early time, 
as explained in App. [Cl positive and negative u wave 
packets follow very different trajectories, whereas for BH 
they essentially follow the same trajectory. 
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